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O ' We report on a spectroscopic study of the spin-wave eigen-modes inside an individual normally 

^*^ ' magnetized two layers circular nano-pillar (Permalloy|Copper|Permalloy) by means of a Magnetic 

I Resonance Force Microscope (MRFM). We demonstrate that the observed spin- wave spectrum crit- 

^fi , ically depends on the method of excitation. While the spatially uniform radio-frequency (RF) 

t^ ' magnetic field excites only the axially symmetric modes having azimuthal index £ = 0, the RF cur- 

rent flowing through the nano-pillar, creating a circular RF Oersted field, excites only the modes 
having azimuthal index £ = +1. Breaking the axial symmetry of the nano-pillar, either by tilting 
^ ' the bias magnetic field or by making the pillar shape elliptical, mixes different ^-index symmetries, 

' which can be excited simultaneously by the RF current. Experimental spectra are compared to 

theoretical prediction using both analytical and numerical calculations. An analysis of the infiuence 
of the static and dynamic dipolar coupling between the nano-pillar magnetic layers on the mode 
spectrum is performed. 
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I. INTRODUCTION 

Technological progress in the fabrication of hybrid 
nanostructures using magnetic metals has allowed the 
emergence of a new science aimed at utilizing spin depen- 
dent effects in the electronic transport properties^. An 
elementary device of spintronics consists of two magnetic 
layers separated by a normal layer. It exhibits the well- 
known giant magneto-resistance (GMR) effect^i^, that is, 
its resistance depends on the relative angle between the 
magnetic layers. Nowadays, this useful property is ex- 
tensively used in magnetic sensors^. The converse ef- 
fect is that a direct current can transfer spin angular 
momentum between two magnetic layers separated by 
either a normal metal or a thin insulating layei^. As a 
result, a spin polarized current leads to a very efficient 
destabilization of the orientation of a magnetic moment*^. 
Practical applications are the possibility to control the 
digital information in magnetic random access memo- 
ries (MRAMs}2,!i2, or to produce high frequency signals 
in spin transfer nano-oscillators (STNOs)ii'i^. 

From an experimental point of view, the precise iden- 
tification of the spin-wave (SW) eigen-modes in hybrid 
magnetic nanostructures remains to be done^^"— . Of 
particular interest is the exact nature of the modes ex- 
cited by a current perpendicular-to-plane in STNOs. 
Here, the identification of the associated symmetry be- 
hind each mode is essential. It gives a fundamental 
insight about their selection rules and about the mu- 
tual coupling mechanisms that might exist intra or in- 
ter STNOs. It also determines the optimum strategy to 
couple to the auto-oscillating mode observed when the 



spin transfer torque compensates the damping, a vital 
knowledge to achieve phase synchronization in arrays of 
nano-pillar&iS,. These SW modes also have a fundamen- 
tal influence on the high frequency properties of these de- 
vices and in particular on the noise of magneto-resistive 
sensors22i2i. 

A natural mean to probe SW modes in hybrid nanos- 
tructures is to use their magneto-resistance properties. 
For instance, thermal SW can be directly detected in the 
noise spectrum of tunneling magneto-resistance (TMR) 
devices owing to their large TMR ratio^i^. It is also pos- 
sible to use spin torque driven ferromagnetic resonance 
(ST-FMR)^"^. In this approach, an RF current flow- 
ing through the magneto-resistive device is used to excite 
the precession of magnetization and to detect it through 
a rectification effect. Direct excitation of SW modes by 
the RF field generated by micro-antennas and their de- 
tection through dc rectification^! or high-frequency GMR 
measurements^ has also been reported in spin- valve sen- 
sors. In all these experiments, the static magnetiza- 
tions in the spin-valve have to be misaligned in order 
for the magnetization precession to produce a finite volt- 
age. Because highly symmetric magnetization trajecto- 
ries do not produce any variation of resistance with time 
in some cases, a third magnetic layer playing the role of 
an analyzer can be introduced^. In ST-FMR, the non- 
collinearity of the magnetizations is also required for the 
RF spin transfer excitation not to vanisb^^i^. More- 
over, the latter was never directly compared to standard 
FMR, where a uniform RF magnetic field is used to ex- 
cite SW modes. Thus, although the voltage detection of 
SW eigen-modes in hybrid nanostructures is elegant, one 



should keep in mind that some of them might be hidden 
due to symmetry reasons. 

Here, we propose an independent metliod of detect- 
ing the magnetic resonance inside a spin- valve nanostruc- 
ture. We shall use a Magnetic Resonance Force Micro- 
scope (MRFM)^"— . A first decisive advantage of the 
MRFM technique is that the detection scheme does not 
rely on the SW spatial symmetry because it measures 
the change in the longitudinal component of the magne- 
tization. Like a bolometric detection, mechanical based 
FMR detects all the excited SW modes, independently 
of their phase^*^. A second decisive advantage is that 
MRFM is a very sensitive technique that can measure the 
magnetization dynamics in nanostructures buried under 
metallic electrodes^"—. Indeed, the probe is a magnetic 
particle attached at the end of a soft cantilever and is 
coupled to the sample through the dipolar interaction. 

In our roadmap to characterize the nature of the auto- 
oscillation modes in STNOs, we report in this work on 
a comprehensive identification of the SW cigcn-modes 
in the simplest possible geometry: the normally mag- 
netized circular spin-valve nano-pillar. This configura- 
tion is obtained by saturating the device with a large 
external magnetic field oriented perpendicular to the lay- 
ers. Thanks to the preserved axial symmetry, a simpli- 
fied spectroscopic signature of the different SW eigen- 
niodcs is expected. This identification is achieved exper- 
imentally from a comparative spectroscopic study of the 
SW eigen-modes excited either by an RF current flow- 
ing perpendicularly through the nano-pillar, as used in 
ST-FMR, or by a homogeneous RF in-plane magnetic 
field, as used in conventional FMR. It shall be developed 
as follows. In section II, we present the MRFM setup 
and the experimental protocol used to perform SW spec- 
troscopy in a spin- valve. We show that the SW spectrum 
excited by a homogeneous RF magnetic field is distinct 
from the SW spectrum excited by an RF current fiowing 
through the nano-pillar. In section III, we perform unam- 
biguous assignment of the resonance peaks to the differ- 
ent layers by experimental means. We determine which 
layer contributes mostly to the observed resonant signals 
by adding a direct current through the nano-pillar, that 
produces opposite spin transfer torques on each magnetic 
layer. In section IV, we analyze the spectra by theoret- 
ical means using both a two-dimensional analytical for- 
malism and a three-dimensional micromagnetic simula- 
tion package, SpinFlow 3D. By careful comparison of the 
measured spectra to the calculations, the nature of the 
SW dynamics in the system is identified and the selec- 
tion rules for SW spectroscopy in perpendicularly mag- 
netized spin-valve nanostructures are established. This 
result is completed in section V by a study of the infiu- 
ence of symmetry breaking on the selection rules. This 
is obtained experimentally by introducing a tilt angle of 
the applied magnetic field, and in simulations by chang- 
ing the shape of the nano-pillar. In the conclusion, we 
emphasize the importance of this work for phase syn- 
chronization of STNOs. The paper is arranged in such 
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FIG. 1. (Color onlme) Schematic representation of the exper- 
imental setup used for this comparative spin-wave spectro- 
scopic study. The magnetic sample is a circular nano-pillar 
comprising a thin Pja and a thick Py^ magnetic layers sep- 
arated by a Cu spacer. It is saturated by a large magnetic 
field Jicxt applied along its normal axis. A cantilever with 
a magnetic sphere attached at its tip monitors the magne- 
tization dynamics inside the buried structure. The inset is 
a microscopy image (top view) of the two independent ex- 
citation circuits: in red the circuit allowing the injection of 
an RF current perpendicular-to-plane through the nano-pillar 
(irf, red arrow); in blue the circuit allowing the generation of 
an RF in-plane magnetic field (/irf, blue arrow). The nano- 
pillar is at the center of the yellow cross-hair. The main figure 
is a section along the A- A direction. 



a fashion so as to present the main results in the body 
of the text. A comprehensive appendix has been put at 
the end of the paper, where the details of the introduced 
material are developed. 



II. FERROMAGNETIC RESONANCE FORCE 
SPECTROSCOPY 

This section starts with a description of the nano-pillar 
sample, followed by a description of the MRFM instru- 
ment used for this spectroscopic study. Then, we com- 
pare the experimental SW spectra excited by an RF cur- 
rent flowing perpendicularly through the nano-pillar, as 
used in ST-FMR, and by a uniform RF magnetic field 
applied parallel to the layers, as used in standard FMR. 



A. The lithographically patterned nanostructure 

The spin-valve structure used in this study is a stan- 
dard Permalloy (Ni8oFc2o=Py) bi-layer structure sand- 
wiching a 10 nm copper (Cu) spacer: the thicknesses of 
the thin Pya and the thick Py^ layers are respectively 
tf, = 4 nm and if, = 15 nm. Special care has been put 
in the design of the microwave circuit around the nano- 
pillar. The inset of FIG. [1] shows a scanning electron 
microscopy top view of this circuit. The nano-pillar is 
located at the center of the cross-hair, in the middle of 



a highly symmetric pattern designed to minimize cross- 
talk effects between both RF circuits shown in blue and 
red, which provide two independent excitation means. 

The nano-pillar is patterned by standard e-beam 
lithography and ion-milling techniques from the extended 
film, (Cu60 I PyaS | CulO | Pya4 | Au25) with thick- 
nesses expressed in nm, to a nano-pillar of nominal ra- 
dius 100 nm. A precise control allows to stop the etch- 
ing process exactly at the bottom Cu layer, which is 
subsequently used as the bottom contact electrode. A 
planarization process of a polymerized resist by reac- 
tive ion etching enables to uncover the top of the nano- 
pillar and to establish the top contact electrode. The top 
and bottom contact electrodes are shown in red tone in 
FIG. [TJ These pads are impedance matched to allow for 
high frequency characterization by injecting an RF cur- 
rent i,f through the device. The bottom Cu electrode is 
grounded and the top Au electrode is wire bounded to the 
central pin of a microwave cable. Hereafter, spectra as- 
sociated to SW excitations by this part of the microwave 
circuit will be displayed in red tone. The nano-pillar is 
also connected through a bias-T to a dc current source 
and to a voltmeter through the same contact electrodes, 
which can be used for standard current perpendicular 
to the plane (CPP-GMR) transport measurements^. In 
our circuit, a positive current corresponds to a flow of 
electrons from the Py^ thick layer to the Pya thin layer 
and stabilizes the parallel configuration due to the spin 
transfer effect^i^. The studies presented below will be 
limited to a dc current up to the threshold current for 
auto-oscillations in the thin layer. 

The originality of our design is the addition of an in- 
dependent top microwave antenna, whose purpose is to 
produce an in-plane RF magnetic field h^f at the nano- 
pillar location. In FIG. [T] this part of the microwave 
circuit is shown in blue tone. The broadband strip-line 
antenna consists of a 300 nm thick Au layer evaporated 
on top of a polymer layer that provides electrical isolation 
from the rest of the structure. The width of the antenna 
constriction situated above the nano-pillar is 10 /im. In- 
jecting a microwave current from a synthesizer inside the 
top antenna produces a homogeneous in-plane linearly 
polarized microwave magnetic field, oriented perpendic- 
ular to the stripe direction. Hereafter, spectra associated 
to SW excitations by this part of the microwave circuit 
will be displayed in blue tone. 



B. Mechanical-FMR 

The nano-fabricated sample is then mounted inside a 
Magnetic Resonance Force Microscope (MRFM), here- 
after named mechanical-FMR^. The whole apparatus 
is placed inside a vacuum chamber (10^^ mbar) oper- 
ated at room temperature. The external magnetic field 
produced by an electromagnet is oriented out-of-plane, 
i.e., along the nano-pillar axis z. The mcchanical-FMR 
setup allows for a precise control, within 0.2°, of the po- 



lar angle between the applied field and z. In our study, 
the strength of the applied magnetic field shall exceed 
the saturation field (rj 8 kOe), so that the nano-pillar is 
studied in the saturated regime. 

The mechanical detector is an ultra-soft cantilever, 
an Olympus Bio-Lever having a spring constant k « 
5 mN/m, with a 800 nm diameter sphere of soft amor- 
phous Fe (with 3% Si) glued to its apex. Standard piezo 
displacement techniques allow for positioning the mag- 
netic spherical probe precisely above the center of the 
nano-pillar, so as to retain the axial symmetry. This 
is obtained when the dipolar interaction between the 
sample and the probe is maximal, by minimizing the 
cantilever resonance frequency, which is continuously 
monitored^. 

The mechanical sensor is insensitive to the rapid oscil- 
lations of the transverse component in the sample, which 
occur at the Larmor precession frequency, i.e., several or- 
ders of magnitude faster than its mechanical resonances. 
The dipolar force on the cantilever probe is thus propor- 
tional to the static component of the magnetization inside 
the sample. For our normally magnetized sample, this 
longitudinal component reduces to Mz. We emphasize 
that for a bi-layer system, the force signal integrates the 
contribution of both layers. Moreover, the local Mz{r) 
in the two magnetic layers is weighted by the distance 
dependence of the dipolar coupling to the center of the 
sphere. In our case though, where the separation be- 
tween the sphere and the sample is much larger than the 
sample dimensions, one can neglect this weighting and 
the measured quantity simplifies to the spatial average: 



{Mz)^^ J^Mz{r)d\, 



(1) 



where the chevron brackets stand for the spatial average 
over the volume of the magnetic body. 

The mechanical-FMR spectroscopy presented below 
consists in recording by optical means the vibration am- 
plitude of the cantilever either as a function of the out- 
of-plane magnetic field i/cxt at a fixed microwave exci- 
tation frequency /gx, or as a function of the excitation 
frequency / at a fixed magnetic field i?fix. This type 
of spectroscopy is called cw, for continuous wave, as it 
is monitoring the magnetization dynamics in the sample 
under a forced regime. A source modulation is applied 
on the cw excitation. It consists in a cyclic absorption 
sequence, where the microwave power is switched on and 
off at the cantilever resonance frequency, /c ~ 11.85 kHz. 
The signal is thus proportional to (AM^), where A rep- 
resents the difference from the thermal equilibrium state. 
The source modulation enhances the signal, recorded by 
a lock-in detection, by the quality factor Q « 2000 of 
the mechanical oscillator. The force sensitivity of our 
mechanical-FMR setup is better than 1 fN, correspond- 
ing to less than 10'^ Bohr magnetons in a bandwidth of 
one second^. We note that this modulation technique 
does not affect the line shape in the linear regime, be- 
cause the period of modulation l//c is very large com- 



pared to the relaxation times of the studied ferromag- 
netic system^i^. Moreover, we emphasize that since 
the mechanical-FMR signal originates from the cyclic 
diminution of the spatially averaged magnetization inside 
the whole nano-pillar synchronous with the absorption of 
the microwave power, it detects all possible SW modes 
without discrimination22ii£. 

Finally, we mention that the stray field produced by 
the magnetic sphere attached on the cantilever does af- 
fect the detected SW spectra. In our setup, the separa- 
tion between the center of the spherical probe and the 
nano-pillar is set to 1.3 //m (see FIG. [T|), which is a large 
distance considering the lateral size of the sample. At 
such distance, the coupling between the sample and the 
probe is weak^^ as it does not affect the profiles of the in- 
trinsic SW modes in the sample. This is in contrast with 
the strong coupling regime, where the stray field of the 
magnetic probe can be used to localize SW modes below 
the MRFM tip^. For our mechanical SW spectrometer, 
the perturbation of the magnetic sphere reduces to a uni- 
form translation of all the peak positions^ by -190 Oe 
(see section |IIIB[) . In the following, all the SW spectra 
are recorded with the magnetic sphere at the same exact 
position above the nano-pillar. 



C. RF magnetic field vs. RF current excitations 

The comparative spectroscopic study performed by 
mechanical-FMR at /fix = 8.1 GHz on the normally mag- 
netized spin- valve nano-pillar is presented in FIG. [H In 
these experiments, there is no dc current flowing through 
the device, and the spectra are obtained in the small ex- 
citation regime (precession angles less than 5°, see ap- 
pendix IB ip . The upper panel (a) shows the SW spec- 
trum excited by a uniform RF magnetic field applied in 
the plane of the layers, while the lower panel (b) dis- 
plays the SW spectrum excited by an RF current flowing 
perpendicularly through the magnetic layers. The strik- 
ing result is that these two spectra are different: none of 
the SW modes excited by the homogeneous RF field is 
present in the spectrum excited by the RF current flow- 
ing through the nano-pillar, and vice versa. 

Let us first focus on FIG. Oi, where the obtained ab- 
sorption spectrum corresponds to the so-called standard 
FMR spectrum. Here, the output power of the microwave 
synthesizer at 8.1 GHz is set to +3 dBm, which corre- 
sponds to an amplitude of the uniform linearly polarized 
RF magnetic field h^f ~ 2.1 Oe produced by the antenna 
(see appendix IB 1|) . In this standard FMR spectrum, 
only SW modes with non-vanishing spatial average can 
couple to the homogeneous RF field excitation. In field- 
sweep spectroscopy, the lowest energy mode occurs at 
the largest magnetic field. So, the highest field peak at 
Hq = 10.69 kOe should be ascribed to the uniform mode. 
Since this peak is also the largest of the spectrum, it 
corresponds to the precession of a large volume in the 
nano-pillar, i.e., the thick layer must dominate in the 
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FIG. 2. (Color online) Comparative spectroscopic study per- 
formed by mechanical-FMR at /fix = 8.1 GHz, demonstrating 
that distinct SW spectra are excited by a uniform in-plane 
RF magnetic field (a) and by an RF current flowing perpen- 
dicularly through the layers (b). The positions of the peaks 
are reported in Table HIl 



dynamics. In mechanical-FMR, a quantitative measure- 
ment of the longitudinal magnetization is obtained^i^ 
(see appendix IB ip . The amplitude of the peak at iJ® 
corresponds to 47r(AA'/^) ~ 14 G, which represents a pre- 
cession angle {9) ^ 3.1°. This sharp peak is followed by a 
broader peak with at least two maxima at iJ® = 9.65 kOe 
and H® = 9.51 kOe, and at lower field, by a smaller reso- 
nance around iJ® = 8.64 kOe. Among these other peaks, 
there is the uniform mode dominated by the thin layer, 
which has to be identified and distinguished from higher 
radial index SW modes. 

Let us now turn to FIG.[2)d, corresponding to the spec- 
troscopic response to an RF current of same frequency 
8.1 GHz flowing perpendicularly through the nano-pillar. 
Here, the output power of the microwave synthesizer is 
-22 dBm, which corresponds to an rms amplitude of 
the RF current i^f ^ 170 fiA (see appendix IB 2p . The 
SW spectrum is acquired under the exact same condi- 
tions as for standard FMR, i.e., the spherical magnetic 
probe of the mechanical-FMR detection is kept at the 
same location above the sample. The striking result is 
that the position of the peaks in FIGS. Oi and [5}d do 
not coincide. More precisely there seems to be a trans- 
lational correspondence between the two spectra, which 
are shifted in field by about 0.5 kOe from each other. 
The lowest energy mode in the RF current spectrum oc- 
curs at Ho = 10.22 kOc. This is again the most intense 



peak, suggesting that the thick layer contributes to it, 
and 4:Tt{AMz) - 26 G, which represents a precession an- 
gle (6) ^ 4.2°. This main resonance line is also split in 
two peaks, with a smaller resonance in the low field wing 
of the main peak, about 100 Oe away. At lower field, 
two distinct peaks appear at H@ = 9.17 kOe and H@ = 
9.07 kOe and another peak is visible at Hq = 8.22 kOe. 

The fact that the two spectra of FIGS. [Hi and[5]D arc 
distinct implies that they have a different origin. It will 
be shown in the theoretical section IIVA3I that the RF 
field and the RF current excitations probe two different 
azimuthal symmetries i. Namely, only £ = modes are 
excited by the uniform RF magnetic field, whereas only 
i = +1 modes are excited by the orthoradial RF Oersted 
field associated to the RF currenii^. The mutually exclu- 
sive nature of the responses to the uniform and orthora- 
dial symmetry excitations is a property of the preserved 
axial symmetry, where the azimuthal index ^ is a good 
quantum number, i.e., different £- index modes are not 
mixed and can be excited separately (see section llV A2I) . 



III. EXPERIMENTAL ANALYSIS 

In this section, we first look at the effect of a continuous 
current flowing through the nano-pillar on the SW spec- 
tra in order to determine which layer contributes mostly 
to the resonant signals observed in FIG. [5] Due to the 
asymmetry of the spin transfer torque in each magnetic 
layer, the different SW modes are influenced differently 
depending on the layer in which the precession is the 
largest. Then, we briefly mention experiments, where 
spectroscopy is performed by monitoring the dc voltage 
produced by the magnetization precession in the hybrid 
nanostructure, and compared to mcchanical-FMR. Fi- 
nally, the analysis of the frequency-field dispersion rela- 
tion and of the linewidth of the resonance peaks enables 
to extract the gyromagnetic ratio and the damping pa- 
rameters in the thick and thin layers. 
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FIG. 3. (Color online) Evolution of the SW spectra measured 
at /fix = 8.1 GHz by mechanical-FMR for different values of 
the continuous current /dc flowing through the nano-pillar. 
The panel (a) corresponds to excitation by a uniform RF mag- 
netic field and the panel (b) to excitation by an RF current 
through the sample. 



A. Direct bias current 

To gain further insight about the peak indexation, we 
have measured the spectral evolution produced on the 
SW spectra of FIG. [D when a finite dc current /dc + 
is injected in the nano-pillar. We recall that for our 
sign convention, a positive dc current stabilizes the thin 
layer and destabilizes the thick one due to the spin trans- 
fer torque, and vice versa^. The results obtained by 
mechanical-FMR are reported in FIG. [3l 

Let us first concentrate on FIG. [3^, in which the ex- 
citation that probes the different SW modes is the same 
as in FIG. [2^, i.e., a uniform RF magnetic field. Two 
main features can be observed in the evolution of the 
SW spectra as /dc is varied. First, the amplitude of the 
peak at //® smoothly increases with the positive current 
and smoothly decreases with the negative current. At 



the same time, the peak at i/®, which is about five times 
smaller than the peak at //® when /dc = mA, almost 
disappears for positive current and strongly increases at 
negative current, until it becomes larger than the other 
peaks when /dc = -4 mA. These two features are consis- 
tent with the effect of spin transfer if we ascribe the peak 
at //{J) to the uniform mode of mostly the thick layer and 
the peak at H@ to the one of mostly the thin layer. More 
precisely, it is expected that in the sub-critical regime 
(|/dc| < ^th, where /th is the threshold current for auto- 
oscillations, /th < for the thin layer and /th > for the 
thick layer), the damping scales as a(l-/dc/^th)^'^ (see 
appendix lA 1[) . where a is the Gilbert damping parame- 
ter. It means that the linewidth of a resonance peak that 
is favored by spin transfer should decrease as the current 
gets closer to /th, and that its amplitude, which scales as 



the inverse linewidth, should increase. 



B. Voltage-FMR 



Although the effect on the peak amplitude noted above 
is clear in FIG.|3^, it is not on the linewidth. The reason 
is that in this experiment, the strength of the driving RF 
magnetic field is kept constant to h^f =2.1 Oe. As a re- 
sult, the shape of the growing peaks in FIG. [5^ becomes 
more asymmetric, which is a signature that the preces- 
sion amplitude driven by the RF field is strong enough 
to change the internal field by an amount of the order 
of the linewidth. This leads to some foldover of the res- 
onance line^i^, a non-linear effect for which details are 
given in the appendix IB II In other words, the distortion 
of the line shape as the peak amplitude increases pre- 
vents to see the diminution of its linewidth^. It would 
be necessary to decrease the excitation amplitude as the 
threshold current is approached^ so as to maintain the 
peak amplitude in the linear regime in order to reveal it. 

The opposite signs of the spin transfer torques which 
influence the dynamics in the thin and thick layers are 
thus clearly seen in FIG. [3^. Their relative strengths can 
also be determined, as the amplitude of the peak at H@ 
grows much faster with negative current than the one of 
the peak at iJ® with positive current. This is because 
the efficiency of the spin transfer torque is inversely pro- 
portional to the thickness of the layeri^. Whereas the 
precession angle in the thick layer does not vary much 
with /dc (from r; 2.5° at -4 mA to Ri 3.5° at +4 mA), 
the precession angle that can be deduced from (Ail/^) in 
the thin layer grows from almost zero at /dc = +4 mA to 
more than 6° at /dc = -4 mA. Moreover, the peak posi- 
tion H@ shifts clearly towards lower field as the negative 
current is increased. This is due to the onset of spin 
transfer driven auto-oscillations in the thin layer, which 
occurs at a threshold current Ith ^ -4 mA and produces 
this non-linear shiflp^. We note, that such a value for 
the threshold current in the thin layer can be found from 
Slonczewski's model (see appendix lAip . 

Let us now briefiy discuss FIG. ^jp, which shows the 
dependence on /dc of the mcchanical-FMR spectra ex- 
cited by an RF current excitation. A similar dependence 
on /dc of the resonance peaks in translational correspon- 
dence with FIG. [5Ji is observed. Again, a clear asymme- 
try is revealed depending on the polarity of /dc and on 
the SW modes. The double peak at //© is favored by 
positive currents, hence it should be ascribed to mostly 
the thick layer precessing, while the double peak at //© 
is strongly favored by negative currents, hence it should 
be ascribed to mostly the thin layer precessing. More- 
over, a careful inspection shows that the peak //©, which 
looks single at /dc = mA, is actually at least double. 
We will explain this splitting of higher harmonics modes 
in section IVBI 

To summarize, the passage of a dc current through the 
nano-pillar enables to determine which layer mostly con- 
tributes to the observed SW modes, owing to the asym- 
metry of the spin transfer effect. 



Our experimental setup also allows to monitor the dc 
voltage produced across the nano-pillar by the precession 
of the magnetization in the bi-layer structure. A lock- 
in detection is used to measure the difference of volt- 
age across the nano-pillar when the RF is on and off: 
Vdc = K)n - VoS- This can be done simultaneously to the 
acquisition of the mechanical-FMR signal, in the exact 
same conditions (see FIG. [1]). Since the presentation of 
the experimental results requires a specific discussion, the 
details as well as the graphs will be published elsewhere. 
Here, we shall only reveal the three main features that 
can be noticed in the voltage-FMR spectra. 

First, even at /dc = 0, dc voltage peaks are pro- 
duced across the nano-pillar at the same positions as 
the mechanical-FMR peaks observed in FIG. [H with a 
difference of potential that lies in the 10 nV range for 
the precession angles excited here. It is ascribed to 
spin pumping and accumulation in the spin-valve hybrid 
structure^^^. Second, these voltage resonance peaks 
arc signed, namely, the SW modes favored at /dc < in 
FIG. [3K (for which the thin layer is dominating) produce 
a positive voltage peak, whereas those favored at /dc > 
(thick layer dominating) produce a negative voltage peak. 
This difference between the thick and thin layer contri- 
butions is ascribed to the asymmetry of the spin accu- 
mulation in the multi-layer stack^. Third, the relative 
amplitudes of the voltage-FMR peaks are different from 
the mechanical-FMR ones. For instance, the voltage- 
FMR peak of the thin layer at H@ is slightly larger than 
the peak at //g) of the thick layer (and it has an opposite 
sign). This illustrates an important difference between 
the two detection schemes. While mechanical-FMR mea- 
sures a quantity proportional to the precessing volume, 
(AMz), the voltage-FMR measures an interfacial effect. 
Therefore, when the same precession angle is excited in 
both layers, the voltage-FMR signal associated to each 
layer is approximately the same, whereas the mechanical- 
FMR signal from the thin layer is roughly four times 
smaller than the one from the thick layer, due to their 
relative thicknesses. 

Finally, we mention that voltage-FMR spectroscopy 
can also record the intrinsic FMR spectrum of the nano- 
pillar, i.e., in the absence of the spherical MRFM probe 
above it. This enables to check that the only effect in- 
troduced by the probe in mechanical-FMR is an overall 
shift of the SW modes spectra to lower field without any 
other distortion, and to quantify this shift, found to be 
-190 Gg^^. 



C. Gyromagnetic ratio 

A precise orientation of the applied magnetic field 
iJcxt along the normal z of the sample (polar angle 
9h = (£,i?oxt) = 0) enables a direct determination of 
the modulus 7 of the gyromagnetic ratio^. By following 



the frequency-field dispersion relation of the resonance 
peaks at Hq and at if® (from 4.5 GHz to 8.1 GHz and 
from 6.2 GHz to 11 GHz, respectively) in our nano-pillar, 
it is found that 7 = 1.87 x 10^ rad.s~^.G"^ is identical in 
the thick and thin layers. Moreover, the value of 7 mea- 
sured in the nano-pillar is the same as in the extended 
reference film (see appendix IB 31 and Table |T| , confirm- 
ing that the applied field is sufficient to saturate the two 
magnetic layers and is precisely oriented along z. 

The same result is obtained by following the evolu- 
tion of the frequency-field dispersion relation presented 
in FIG. m Here, we take advantage of the broadband 
design of the electrodes which connect the nano-pillar 
to measure the FMR spectrum at fixed bias magnetic 
field, Hfix = 10 kOe, by sweeping the frequency of the 
RF current through it. The data are plotted accord- 
ing to the frequency scale above FIG. |4^. At con- 
stant magnetic configuration (above the saturation field, 
i.e., > 8 kOe), this frequency scale is in correspondence 
with field-sweep experiments performed at fixed RF fre- 
quency /fix = 8.1 GHz through the affine transformation 
Hcxt - Hfix = 27r(/ - /fix)/7, as seen from the field scale 
below FIG. |31d. This is a direct experimental check of 
the equivalence between frequency and field sweep ex- 
periments in the normally saturated state. 



D. Damping parameters 

From the FMR data presented above, we can also di- 
rectly extract the damping parameters in each Permalloy 
layer. Indeed, in field-sweep spectroscopy in the normal 
orientation {9h = 0), the full width at half-maximum 
(FWHM) AH of a resonance line is proportional to the 
excitation frequency aj/(27r) through the Gilbert constant 
a: AH = 2a{uj/j) (see appendix I A ip . 

The linewidth of the peak at Hq associated to mainly 
the thick layer in FIG. [5^ is equal to AHq = 48 Oe, 
which corresponds to a damping a® = 0.88 x 10^^. From 
the same mechanical-FMR spectrum, the linewidth of 
the peak at fl@, associated to mainly the thin layer, can- 
not be easily extracted due to the proximity of the peak 
at iJ®. Owing to the interfacial origin of the voltage- 
FMR signal, the peak at H@ is more distinguishable in 
the spectrum of the voltage-FMR (not shown), and its 
linewidth, AH@ = 70 Oe, can be fitted. It corresponds to 
a damping a® = 1.29 x 10"^. 

The linewidths of the modes at Ho and iJ© can also be 
fitted and give similar results for the damping associated 
to each layer. In the case of the RF current excitation, 
a frequency-sweep spectrum can be acquired at a fixed 
bias magnetic field H^x (see FIG. |3|). In that case, the 
damping constant is simply obtained by a = A//(2/), 
where A/ is the width of the line centered at /. At 
iJflx = 10 kOe, /o = 7.37 GHz and A/o = 0.12 GHz, 
which yield a© = 0.81 x 10"^ and /© = 10.92 GHz and 
A/© = 0.33 GHz, which yield a© = 1.5 x 10"^ 

In summary, we retain the following values for the 
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FIG. 4. (Color online) Frequency-field dispersion relation: the 
top spectrum (a) is measured at fixed bias field Ha^ = 10 kOe 
by sweeping the frequency of the RF current irt through the 
nano-pillar. The bottom spectrum (b) is the same as in 
FIG. ^jp, and is obtained by sweeping the magnetic field at 
fixed frequency /ax = 8.1 GHz of Vf. The top and bottom 
scales are in correspondence through the affine transforma- 
tion Hoxt - -fffix = 27r(/ - /flx)/7- 



damping parameters in respectively the thin and the 
thick layers: aa = (1.4 ± 0.2) x 10"^ and ab = (0.85 ± 
0.1) X 10"^. We have reported them, together with 7. in 
Table HI 

These two values are in line with the ones obtained on 
the reference film, which have also been reported in Ta- 
ble |T1 Still, we observe that the linewidths in the nanos- 
tructure are systematically lower than the ones measured 
on the reference film. This is a constant characteris- 
tic that we associate to the confined geometry, which 
lifts most of the degeneracy (well separated SW modes) 
and thus strongly reduces the inhomogeneous part of the 
linewidth observed in the infinite laye r^^'^^ . Rather, the 
inhomogeneities associated to the magnetic layersi^ or 
to the confinement geometry will lead to some mode 
splitting in the nanostructure (see section IVB[) . We 
have checked that the inhomogeneous contribution to the 
linewidth in the nano-pillar is weak, by following the de- 
pendence of the measured AH as a function of frequency. 
In fact, the increase of AH@ from 70 Oe at 8.1 GHz to 
105 Oe at 11 GHz is purely homogeneous. 

Finally, the finding that the damping is larger in the 
thin layer than in the thick layer is ascribed to the ad- 
jacent metallic layers^. In fact, non-local effects such 



TABLE I. Magnetic parameters of the thin Pya and thick 
Pyi, layers measured by cavity-FMR on the reference film 
(top row) and by mechanical-FMR in the nano-pillar (bot- 
tom row). 
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as the spin pumping effec t^ '^^ and the spin diffusion in 
the adjacent normal layers by the conduction electrons 
yield an interfacial increase of the magnetic damping^, 
stronger in the case of thin layers. 



IV. THEORETICAL ANALYSIS 

In this section, we first review a general formalism al- 
lowing the calculation of the discrete spectrum associ- 
ated with SW propagation inside a confined body of ar- 
bitrary magnetic configuration. It is shown that in the 
two-dimensional (2D) axially symmetric case, different 
£-index modes can be excited separately, as found ex- 
perimentally in section III CI The classification of the 
SW modes in this case is also used to extract the pa- 
rameters of each magnetic layer from the experimental 
FMR spectra. In a second part, we discuss the influ- 
ence of the dynamic coupling between the magnetic disks, 
where the collective dynamics splits into binding and 
anti-binding modes. It is shown that in our experimen- 
tal case, the dynamic dipolar coupling introduces a weak 
spectral shift, although its influence on the character of 
the SW modes is real. In the last part, a comparison to 
full three-dimensional (3D) micromagnetic simulations is 
performed in order to study in details the collective dy- 
namics in the nano-pillar. 



A. Analytical model 

1. General theory 

Below, we briefly review the general theory of linear 
SW excitations (see appendix lA II for more details). We 
consider an arbitrary equilibrium magnetic configuration, 
where the local magnetization writes MgU, with Ms the 
saturation magnetization and u the unit vector along the 
local equilibrium direction (implicitly dependent on the 
spatial coordinates). The linearization of the local equa- 
tion of motion is obtained by decomposing the instanta- 
neous magnetization vector M{t) into a static and dy- 
namic componenli^ (see FIG. O. We shall use the fol- 




FIG. 5. (Color online) Schematic representation of the mag- 
netization dynamics under continuous RF excitation. In the 
steady state, the torque exerted by the RF perturbation field 
hi (orange arrow) compensates the torque exerted by the 
damping (green), and the local magnetization vector M{t) 
(purple) precesses at the Larmor frequency on a circular orbit 
around the local equilibrium direction (unit vector u). M{t) 
is the vector sum of a small oscillating component Msm and 
a large static component Ms (l - \mf /2), respectively trans- 
verse and parallel to ii. The inset shows the simulated spa- 
tial distribution of u inside the nano-pillar at iJcxt = 10 kOe 
(see section riVC[) . In the white regions, the magnetization 
is aligned along the normal z within 0.05°. In the colored 
regions, u is fiaring (< 5°) in the radial direction (the hue 
indicates the direction of ■u - z according to the color code 
defined in FIG. [6|. 



lowing ansatz: 



Mjt) 



u + m{t) +0(m^). 



(2) 



where the transverse component m{t) is the small dimen- 
sionless deviation (|97i| « 1) of the magnetization from 
the equilibrium direction. In ferromagnets, |7Vf | = M^ is 
a constant of the motion, so that the local orthogonality 
condition -u ■ m = is required. 

Substituting Eq. ^ in the lossless Landau-Lifshitz 
equation Eq. (|A1|) (see appendix I A ![) and keeping only 
the terms linear in m, one obtains the following dynam- 
ical equation for m: 



dm 

~dt 



■u X Sl =(- m, 



(3) 



where here and henceforth, tensor operators are indicated 
by wide hat, the cross product is denoted by x and the 
convolution product is denoted by x-. The self-adjoint 
tensor operator Jl represents the Larmor frequency: 



rj = 7iJ/ + 47r7M,G, 



(4) 



where 7 is the modulus of the gyromagnetic ratio, H 
is the scalar effective magnetic field, / is the identity 
matrix, and G is the linear tensor operator describing 



the magnetic self-interactions. The later is the addi- 
tion of several contributions G'-''-' + G^'^-' + ..., respectively 
the magneto-dipolar interactions, the inhomogeneous ex- 
change, etc... (see appendix IA2p . The effective magnetic 
field H is a, vector aligned along u, whose norm is 



H = ii- Hf)- AttMsU • G * m , 



(5) 



the sum of the ■u-component of Hq^ the total applied 
magnetic field including the stray field of any nearby 
magnetic object (in our case, the adjacent magnetic layer 
in the nano-pillar and the spherical probe), reduced by 
the static self-interactions, which include the depolariza- 
tion magnetic field along u created by the static compo- 
nent of the magnetization. 

SW modes m^ are by definition eigen-solutions of 
Eq. ©: 



■ iuj,,m,, = uxfl * m,. 



(6) 



Here uJu is the SW eigen- frequency and z^ is a set of indices 
to enumerate the different modes. 

The main properties of SW excitations follow from the 
eigen problem Eq. ^ and the fact that the operator S7 
is self-adjoint and real. One can show that the eigen 
solutions obey the closure relation 



i(m^ ■ {u X rriy,)) = NyS^^, 



(7) 



where 5 is the Kronecker delta function and rn stands 
for the complex conjugate of m. Here we have used 
the chevron bracket notation introduced in Eq. ([T]) to 
denote the spatial average. The quantities M^ are real 
normalization constants, which depend on the choice of 
eigen- functions m^ . If the equilibrium magnetization u 
corresponds to a (local) minimum of the energy, then the 
operator Q, is positive-definite. It follows that the "phys- 
ical" modes with lo^ > Q have positive norm Mu > 0. In 
this formalism, the eigen-frequencies ui^ can be calculated 
as 



UJy 



{nil, ■ ^ * rriv) 



(8) 



The importance of this relation is that the frequencies 
LUy calculated using Eq. ([8]) are variationally stable with 
respect to perturbations of the mode profile rriy. Thus, 
injecting some trial vectors inside Eq. ^ allows one to 
get approximate values of uji, with high accuracy^. The 
trial vectors should obey some simple properties: i) they 
should form a complete basis in the space of vector func- 
tions m, ii) be locally orthogonal to u and iii) satisfy 
appropriate boundary conditions at the edges of the mag- 
netic body^. 



2. Normally magnetized disks 

In this part, we shall establish a SW modes basis m^ 
for a normally magnetized disk. A specific feature of the 



considered geometry is its azimuthal symmetry. Mathe- 
matically, this means that the operator m x Jl commutes 
with the operator Rz that describes an infinitesimal ro- 
tation about the z axis, assuming that the boundary con- 
ditions are invariant under such a rotation. 

This particular configuration allows us to classify the 
SW modes according to their behavior under the rota- 
tions in the (a;,y) plane. Namely, SW cigen-modcs are 
also eigen-functions of the operator Rz corresponding to 
a certain integer azimuthal number £: 



dm 

"90 



z X 171= -i{£- l)m. 



(9) 



Here, (f> is the azimuthal angle of the polar coordinate 
system. 

As one can see, Eq. (j9]) determines the vector struc- 
ture of SW modes and their dependence on the angle (p. 
Namely, Eq. ([9]) for a fixed £ has two classes of solutions: 



.(1) 



1 



mY' = -{x + iy)e-"''^i^^/>{p) 



-^^'p.iMh 



(lOa) 



and 



m 



f = i(&-*y)e-(^-2)Vr(p), (10b) 



/,(l-2) 



where the functions ipi, ' (p) describe the dependence of 
the SW mode on the radial coordinate p and have to be 
determined from the dynamical equations of motion. So, 
the azimuthal symmetry allows one to reduce the 2D (p 
and (f>) vector equations to a one-dimensional (p) scalar 
problem. 

Generally speaking, SW eigen-modes are certain lin- 
ear combinations of both possible £-forms Eqs. pUj) . The 
coupling of these two forms is due solely to the inho- 
mogeneous dipolar interaction. In our experimental case 
(lowest energy modes of a relatively thin disk) one can 
completely neglect this coupling^ and consider only the 
right-polarized form Eq. (|10a|) . In the following we will 
drop the superscript (1) in m^ and ip^ . 

Wc shall now find an appropriate set of radial func- 
tions ipe{p) to calculate the SW spectrum using Eq. ([8|). 
Here, we can take advantage of the variational stability 
of Eq. ([5]) and, instead of the exact radial profiles ipe{p) 
(to find them one has to solve integro-differential equa- 
tions), use some reasonable set of functions. Namely, 
it is known that the dipolar interaction in thin disks or 
prisms does not change qualitatively the profile of SW 
modes, but introduces effective pinning at the lateral 
boundaries^. Therefore, we will use radial profiles of the 
form ipi{p) = Ji{ki^nP), where Je{x) is the Bessel func- 
tion and ki^n are SW wave-numbers determined from the 
pinning conditions at the disk boundary p = R. For our 
experimental conditions {ta,tb « -R), the pinning is al- 
most complete, and we shall use k^^n = i^i,nlR, where 
Ki^n is the n-th root of the Bessel function of the ^-th 
order. 
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FIG. 6. (Color online) Color representation of the Bessel spa- 
tial patterns for different values of the azimuthal mode index 
£ (by row) and radial mode index n (by column). The arrows 
are a snapshot of the transverse magnetization m^, labeled 
by the index u = £,n. All arrows are rotating synchronously 
in-plane at the SW eigen- frequency. In our coding scheme, 
the hue indicates the phase </> = arg(Tni,) (or direction) of 
m^, and the brightness the amplitude of \mi,\ . The nodal 
positions {\m,^\ = 0) are marked in white. 



FIG. [S] shows a color representation of the Bessel spa- 
tial patterns for different values of the index ly = £,n. 
We restrict the number of panels to two values of the 
azimuthal mode index, i = 0,+l, with the radial index 
varying between n = 0,1,2. In our color code, the hue 
indicates the phase (or direction) of the transverse com- 
ponent m^, while the brightness indicates the amplitude 
of Irriiyp. The nodal positions are marked in white. A 
node is a location where the transverse component van- 
ishes, i.e., the magnetization vector is aligned along the 
equilibrium axis. This coding scheme provides a distinct 
visualization of the phase and amplitude of the preces- 
sion profiles. The black arrows are a snapshot of the m^ 
vectors in the disk and are all rotating synchronously in- 
plane at the SW eigen-frequency. 

The top left panel shows the z^ = 0,0 (^ = 0, n = 0) 
mode, also called the uniform mode. It usually corre- 
sponds to the lowest energy mode since all the vectors 
are pointing in the same direction at all time. Below is 
the £ = +1, n = mode. It corresponds to SWs that 
are rotating around the disk in the same direction as 
the Larmor precession. The corresponding phase is in 
quadrature between two orthogonal positions and this 
mode has a node at the center of the disk. The variation 
upon the n = 0,1,2 index {i being fixed) shows higher 
order modes with an increasing number of nodal rings. 
Each ring separates regions of opposite phase along the 
radial direction. All these spatial patterns preserve the 
rotation invariance symmetry. 
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FIG. 7. (Color online) Analytically calculated spectra at 
Hflx = 10 kOe using the set of Bessel functions (see FIG. |6]) 
as the trial eigen- vectors. The panel (a) shows the linear re- 
sponse to a uniform excitation field hi = x and the panel (b) 
to an orthoradial excitation field hi = -sin<j)x + cos (jjy. A 
light (dark) color is used to indicate the energy stored Eq. (|12|) 
in the thin Py^ and thick Py^ layers. 



3. Selection rules 

Using the complete set of Bessel functions in Eq. ([S|), 
one can obtain analytically the discrete spectrum of 
eigen-values for both the thin and thick layers. The de- 
tails of the numerical application can be found in ap- 
pendix |X2j The spectral values are displayed in FIG. [7] 
using vertical ticks labeled v = jin, where j = a,b indi- 
cates the precessing layer, and £, n the azimuthal and 
radial mode indices. They are calculated at fixed applied 
field i?fix = 10 kOe and placed on the graphs according 
to the frequency scale below FIG. [7|d, which is in corre- 
spondence with the field scale above FIG. [7^ (see IIIICI 
for the equivalence between field- and frequency-sweep 
experiments). 

The comparison with the experimental data in 
FIGS. [2^ and [2}j shows that the coupling to an exter- 
nal coherent source depends primarily on the ^-indcx. 
Indeed; this index carries the discriminating symmetry 
in SW spectroscopy^. This is because the excitation ef- 
ficiency is proportional to the overlap integral 



K 



{m„-hi) 



(11) 



where hi{r) is the spatial profile of the external cxci- 
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tation field. It can be easily shown that a uniform RF 
magnetic field, hi = h^fX, can only excite i= SW modes. 
We have shown in FIG. [7^ the predicted position of these 
modes with blue tone ticks. Obviously the largest over- 
lap is obtained with the so-called uniform mode (n = 0). 
Higher radial index modes (n t 0) still couple to the uni- 
form excitation but with a strength that decreases as n 
increase o^^'^^ . The £+ normal modes, however, are hid- 
den because they have strictly no overlap with the excita- 
tion. The comparison with the experimental spectrum in 
FIG. [5^ confirms that conventional FMR^ probes only 
partially the possible SW eigen-modes, along the i = 0- 
index value. In contrast, the RF current-created Oersted 
field, hi = hoc{p){-sm(j)x + cos4>y) has an orthoradial 
symmetry and can only excite i = +1 SW modes. We 
have shown in FIG. [TJa the predicted position of these 
modes with red tone ticks. They are in good agreement 
with the resonance positions observed in FIG. [Dd. We 
also note that the i = and i = +1 spectra calculated 
analytically bear similar a/b and n index series as a func- 
tion of energy. This explains why the two spectra in 
FIGS. [Da, and [5]d look in translational correspondence 
with each other. We emphasize that the same transla- 
tional correspondence would have been observed for any 
higher azimuthal order i > 1 index spectra. 

From the coupling to the excitation field expressed by 
Eq. pTj) . one can also calculate the mechanical-FMR 
signal oc (AAf^), proportional to the energy stored in 
the magnetic system^Si^. For an arbitrary pulsation fre- 
quency OJ, 



4tt {AM ■ u) c^ AttMsY, 



1^\K 



7f (c.-c..)2 + r2 



■M,., 



(12) 



where the SW damping rate T^ is given by Eq. ()A8|) in ap- 
pendix |XTJ Eq. (|12p is derived under the approximation 
that the only relevant coefficients in the damping matrix 
are the diagonal terms. It has been used to compute the 
relative peak amplitudes in the analytically calculated 
spectra of FIG. H 



4- Comparison with experiments 



The analytical model outlined in sections IIV A II and 
IIV A2I can be used to analyze the experimental spec- 
tra of FIG. m and to extract some useful parameters 
of the nano-pillar. More details can be found in the 
appendix IA2I along with an approximate expression for 
the SW frequencies in the form of Kittel's traditional 
formula (with renormalized values of the effective self- 
demagnetization fields). This Kittel's formula, derived 
for the ^ = spectrum, should be used to analyze the 
SW spectrum excited by a uniform RF field to yield the 
correct values of the magnetization in our nano-pillar. 
Identifying the experimental peaks at H@ and i?® as the 
lowest energy modes of the thin Pya and thick Py;, layers 
yields their respective magnetizations AnMa = 8.0 x 10^ G 



and An Mb = 9.6x 10^ G, sec Eq. (|A32|) . These values have 
been reported in Table HI together with those measured 
in the reference film (see appendix IB3[) . The magneti- 
zations extracted in the nano-pillar are the same as in 
the extended film. The only small difference concerns 
the magnetization of the thin layer, which is 200 G lower 
in the nanostructure than in the reference film (where 
AirMa = 8.2 X 10^ G). We attribute this to some interdif- 
fusion between Py and Cu or Au at the interfaces of the 
thin layer, which can happen during the etching process 
of the nano-pillar. 

Second, the separation between SW modes crucially 
depends on the lateral confinement in the nano-pillar and 
thus on the precise value of its radius. Experimentally, 
the measured field separation between the two first peaks 
in FIG. Ha, (FIG. Ha), which differ by an additional node 
in the radial direction, is Hq - H® = 1.04 kOe {Ho - 
H@ = 1.05 kOe). Using the nominal radius 100 nm in the 
analytical model predicts that consecutive n-indcx mode 
(n = and n= 1 modes) should be separated by 1.33 kOe, 
which is larger than the observed value. This separation 
drops to 1.05 kOe for a larger disk radius R = 125 nm, 
which we thus refer to as the radius of our nano-pillar. 
This value of R also allows to estimate the shift between 
the £ = and £ = +1 spectra, found to be 530 Oe, in good 
agreement with the experimental value H(X)-Hq = 470 Oe 
observed in FIG. H 



B. Influence of dipolar coupling between different 
layers 

In the treatment above we have neglected the dynamic 
coupling between the two magnetic disks in dipolar inter- 
action. In general, the interaction between two identical 
magnetic layers will lead to the hybridization of the same 
jz-index mode of each layer into two collective modes: the 
acoustic mode, where the layers are precessing in phase, 
and the optical mode, where they are precessing in anti- 
phase. This has been observed in interlayer-exchange- 
coupled thin films^ and in trilayered wires where the 
two magnetic stripes are dipolarly coupled^. In the case 
where the two magnetic layers are not identical (different 
geometry or magnetic parameters), this general picture 
continues to subsist. Although both isolated layers have 
eigen-modes with different cigcn-frcquencies, the collec- 
tive magnetization dynamics still splits in a binding and 
anti-binding state. But here, the precession of magneti- 
zation can be more intense in one of the two layers and 
the spectral shift of the coupled SW modes with respect 
to the isolated SW modes is reduced, as it was observed 
in both the dipolarly-^ and exchange-coupled cascs^. 

Here, we assume that the dominant coupling mecha- 
nism between the Py layers is the magnetic dipolar inter- 
action. We neglect any exchange coupling between the 
magnetic layers mediated through the normal spacer or 
any coupling associated to pure spin currenta" in our all- 
metallic spin- valve structure. To analyze the infiuence of 
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FIG. 8. (Color online) Schematic representation of the cou- 
pled dynamics between two different magnetic disks. Here, 
Ldb, the eigen- frequency of the lowest energy precession mode 
in the thick layer (the thin layer being fixed at equilibrium) 
is smaller than tOa, the one in the thin layer (the thick layer 
being fixed at equilibrium). When the two disks are dynam- 
ically coupled through the dipolar interaction, the binding 
state B corresponds to the two layers oscillating in anti-phase 
at lob, with the precession occurring mostly in the thick layer, 
whereas the anti-binding state A corresponds to the layers os- 
cillating in phase at uja, with the precession mostly in the thin 
layer. This is shown by displaying the dipolar charges and the 
precession profile m{p) in each layer using a light (dark) color 
to represent the contribution of the thin (thick) layer. 



the dipolar coupling between the two magnetic layers, 
one can complement the perturbation theory derived in 
the previous section HV Al and in the appendix lA 11 De- 
noting Cj, the SW amplitudes in j-th disk , one can get 
from Eq. (jAe]) : 



dcg 
dt 



-ILOaCa + ilK.bCb , 



dcb . 

-— = -lUJbCb + Ijh^aCa ■ 

at 



(13a) 
(13b) 



where ujj is the frequency of the j-th disk (j = a, b) with 
account of only the static field of the j'-th disk (j' = b, a) 
{i.e., with M'j fixed at equilibrium, see FIG. [3). The 
cross term hjj' is given by 



AttM 



'jj 



M 



^{rn,-G^^Km,A 



(14) 



Here, G'*^-* represents the magneto-dipolar interaction, 
Mj/ is the saturation magnetization of the j'-th disk and 
the averaging goes over the volume of j-th disk. Thus, 
hjji is the average over the j-th mode of the magnetic 
field created by the magnetization of the j'-th disk. It 
can be shown that the overlap defined in Eq. ffT^ is max- 
imum between mode pairs bearing similar wave-numbers 
in each layer (i.e., the same set of indices i/)^. This is 
the reason why dropping the index v in Eqs. (jl3|) and 
p4[) is a reasonable approximation. 



The anti-binding (A) and binding (B) 
frequencies of Eqs. p^ have the form 



eigen- 



^A,B 



OJa + Wfc 



/(Uja -^bV 



+ n'^, 



where 



r2 = 7 ha.bhb.a 



(15) 



(16) 



In the case when the dipolar coupling is small {Q, « 
\uja - wti|), the eigen- frequencies can be written as (we 
assume uja > iOh) 

LOA = ^a+ , (17) 



iOa -UJb 



Wb = Wf, ■ 



Wa 



■UJb 



(18) 



These equations can be used for quantitative purposes 
when rj/lwf, - LJbl < 0.3 in which case they describe fre- 
quency shift with accuracy better than 10%. Thus, the 
larger of the frequencies (wj,) shifts up by 



n' 



-UJb 



(19) 



while the smaller one (ujb) shifts down by the same 
amount. This effect is summarized in FIG. [8] 

A numerical estimate of the coupling strengths ha^b 
and hb^a between the lowest energy SW modes in each 
disk can be found in appendix IA2I The obtained re- 
sult is very close to the approximate estimation used in 
Refjii, where the spatial structure of the interacting SW 
modes is ignored to calculate the dipolar coupling be- 
tween uniformly precessing disks. For the experimental 
parameters, VLJ2-K c^ 0.5 GHz. This coupling is almost an 
order of magnitude smaller than the frequency splitting 
uja-uJb, caused, mainly, by the difference of effective mag- 
netizations of two disks: ^^^{Mb - Ma) - 2tt • 4.5 GHz. 
As a result, the shift of the resonance frequencies due to 
the dipolar coupling is negligible, Au!/2tt ~ 0.06 GHz. 

Using Eqs. (fT5|) . one can also estimate the level of mode 
hybridization due to the dipolar coupling. For instance, 
at the frequency uja ~ Wa, the ratio between the preces- 
sion amplitudes in the two layers is given by 



\cb/c 



Auj/{jha,b) 



n 



■UJb 



(20) 



For the experimental parameters, fl/{uJa -ujb) ~ 0.1, i.e., 
the precession amplitude in the disk b is about 10 % of 
that in the disk a. Thus, although the dipolar coupling 
induces a small spectral shift (second order in the cou- 
pling parameter, Eq. (|19p ). its influence in the relative 
precession amplitude is significant (first order in the cou- 
pling parameter, Eq. ()20p ). Finally, we point that here 
the dipolar coupling is anti-ferromagnetic, and that the 
binding (lower energy) mode B always corresponds to 
the thick layer mainly precessing, with the thin layer vi- 
brating in anti-phase, and vice- versa for the anti-binding 
(in-phase) mode A (see dipolar charges in FIG.[S]). 
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FIG. 9. (Color online) Panel (a) is the numerically calculated 
spectral response to a uniform excitation field hi oc x, from 
a 3D micromagnetic simulation performed at -ffflx = 10 kOe. 
The peaks are labeled according to their precession profiles 
shown in FIG. 1111 A light (dark) color is used to indicate the 
energy stored in the thin (thick) layer. Panel (b) recalls the 
experimental spectrum measured by mechanical-FMR when 
exciting the nano-pillar by a homogeneous RF magnetic field 
at /fix = 8.1 GHz. 



C. Micromagnetic simulations 

In the analytical formalism presented above, several 
approximations have been made. For instance, we have 
assumed total pinning at the disks boundary for the SW 
modes and no variation of the precession profile along 
the disks thicknesses (2D model), and we have neglected 
the dependence on u of the dynamic dipolar coupling. 
Still, it allows to extract important parameters in our 
nano-pillar, such as its radius and the magnetization in 
both layers. It also describes the influence of the dynamic 
dipolar coupling on the position and collective character 
of the SW modes. 

Instead of developing a more complex analytical for- 
malism, we have performed innovative 3D micromagnetic 
simulations in order to go beyond the approximations 
mentioned above, and to unambiguously identify the SW 
modes observed in our nano-pillar sample. For that pur- 
pose, we have used a combination of micromagnetic sim- 
ulation solvers available as part of SpinFlow 3D, a finite 
element based simulation platform for spintronics devel- 
oped by In Silicio^. The steady state micromagnetic 
solver used to obtain numerical approximations of mi- 
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FIG. 10. (Color online) Panel (a) is the simulated spectral 
response to an orthoradial excitation field hi oc -sin<^£ + 
cos (j) y. Panel (b) recalls the experimental spectrum measured 
by mechanical-FMR for an RF current excitation. 



cromagnetic equilibrium states is based on a weak for- 
mulation and Galerkin type finite element implementa- 
tion of the very efficient projection scheme introduced 
in Rcf.— . A second numerical solver, a micromagnetic 
Eigen solver, has been used for fast calculations of loss- 
less 3D SW eigen-modes. It is based on a finite element 
discretization of the generalized eigen- value problem de- 
fined by the linearized lossless magnetization dynamics 
in the vicinity of an arbitrary pre-computed equilibrium 
state, following an approach very similar to the one in- 
troduced in Ref.>2i. The discrete generalized eigen- value 
problem is solved with an iterative Arnoldi method us- 
ing the ARPACK library^^. In this calculation the full 
complexity of the 3D micromagnetic dynamics of the 
presently considered bilayer system is preserved. The 
solver outputs both the eigen-values by increasing en- 
ergy order and the associated eigen- vectors. Several tens 
of SW eigen-modes can be accurately computed in a mat- 
ter of few minutes of CPU time with a standard desktop 
PC, for magnetic thin film nano-structures with typi- 
cal lateral sizes in the 100 nm range. This is two to 
three orders of magnitude faster compared to the re- 
quired computation time when using more traditional 
approaches for micromagnetic computation of SW eigen- 
modes, which are typically based on the Fourier com- 
ponent analysis of time series generated by the solution 
of the full non-linear Landau-Lifshitz-Gilbert equation^. 
Finally, a quite generic linear response solver, implement- 
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ing among other things the spectral decomposition of the 
MRFM signal as expressed in Eqs. (dH, ^ and ([M]) . 
has been used to compute the MRFM spectra shown 
here. 

To proceed, the nano-pillar is first discretized using 
unstructured meshing algorithms resulting in an aver- 
age mesh size of 3.5 nm. This corresponds to a total 
number of vertices in the vicinity of 5 x 10**. The mag- 
netization vector is interpolated linearly inside each cell 
(tetrahedra) - a valid approximation taking into account 
that the cell sizes are smaller than the exchange length 
Acx - 5 nm in Permalloy. The magnetic parameters in- 
troduced in the code are the ones reported in Table HI 
and the simulation incorporates the perturbing presence 
of the magnetic sphere attached on the cantilever. More- 
over, the 10 nm thick Cu spacer is replaced by vacuum, 
so that the layers are only coupled through the dipolar 
interaction (spin diffusion effects are absent). 

The next step is to calculate the equilibrium config- 
uration in the nano-pillar at i/cxt = ^fix = 10 kOe. 
The external magnetic field is applied exactly along z 
and the spherical probe with a magnetic moment m = 
2 X 10"^° emu is placed on the axial symmetry axis at a 
distance s = 1.3 /im above the upper surface of the nano- 
pillar. The convergence criterion introduced in the code 
is \dMz/Mj\ < 2 ■ 10"^ between iterations. The result 
shown in the inset of FIG. [5] reveals that the equilib- 
rium configuration is almost uniformly saturated along 
z. Still, a small tilt (< 5°) of the magnetization, away 
from z and along the radial direction, is observed at the 
periphery of the thick and thin layers. 

The micromagnetic eigen solver is then used to com- 
pute the lowest eigen- values of the problem as well as the 
associated eigen- vectors. The discrete list of eigen- values 
under 18 GHz is shown as black vertical ticks at the bot- 
tom of FIGS. [9^ andfTOk. The precession patterns of the 
six eigen-vectors corresponding to the six lowest eigen- 
frequencies are shown in FIG. [TTJ The middle and right 
columns show the dynamics m in the thin Py^ and thick 
Pyfc layers, while the precession profiles along the median 
direction are shown on the left in light and dark colors, 
respectively. The resonance peaks are labeled according 
to the SW modes precession profiles and the eigen- values 
of the simulated peaks are reported in Table [Til 

From the eigen-vectors spatial patterns, one can com- 
pute their coupling (Eq. (|TT|) ) to a uniform RF field hi = 
hrix and, with Eq. p^. the mechanical-FMR spectrum 
(FIG. ^). The same procedure is repeated for the RF 
current-induced Oersted field hi ex {^.f (^- sin cf) x + cos (j) y) 
excitation fFIG.fTUb). Since the code gives access to the 
contribution of each layer, a light (dark) tone is used 
to indicate the vibration amplitude in the thin (thick) 
layer in the two figures. For comparison, the mechanical- 
FMR spectra of FIGS. [2ti andU^ have been reported in 
FIGS, int and [TUb . respectively. We have applied the 
same conversion between the frequency (top) and field 
(bottom) scales as discussed in section UlI CI 

In FIG. ink, the largest peak in the simulation occurs at 



TABLE II. Comparative table of the resonance values for the 
SW modes, arranged in order of increasing energy. On the left 
are the consecutive peak locations measured experimentally. 
Experiments are performed at /ax = 8.1 GHz (FIG. ^ or 
Hflx = 10 kOe (FIG. Ilk). On the right are the simulated 
eigen-frequencies / at Ha^ = 10 kOe. The conversion to field 
value -ffoxt is obtained through ffcxt - H^^ = 27r(/ - /flx)/7. 



Exp. 


/ (GHz) 
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the same field as the experimental peak at iJ®. This low- 
est energy mode corresponds to the most uniform mode 
with the largest wave-vector and no node along the ra- 
dial direction, thus it has the index n = 0. It has uniform 
phase along the azimuthal direction, which is the char- 
acter of the ^ = index. For this mode, the thick layer is 
mainly precessing, with the thin layer oscillating in anti- 
phase (binding index B), as can be seen from its spatial 
profile in FIG. HU The same analysis can be made for 
the second peak, labeled Bqi, which occurs close to the 
peak at H®. It also corresponds to a resonance mainly 
of the thick layer, and its color representation shows that 
this is the first radial harmonic (n = 1), with one line of 
nodes in the radial direction. Again, the thin layer is os- 
cillating in anti-phase, with the same radial index n = 1, 
as clearly shown by the mode profile along the median 
direction. The third peak is labeled Aqq and is located 
close to the experimental peak at iJ®. It corresponds this 
time to a uniform (n = 0) precession mainly located in the 
thin layer, in agreement with the experimental analysis 
presented in section llll Al In this mode, the thick layer is 
also vibrating in phase with the thin layer (anti-binding 
index A). 

We can also look at the relative amplitudes of preces- 
sion in the two disks to quantify the dynamic coupling 
between the disks. From the profiles shown in FIG. [TTl 
one can infer that for the fundamental mode Bqq, the am- 
plitude of precession is distributed with a ratio of about 
3:1 between the thick (75%) and the thin layer (25%). 
For the mode ^oo, the ratio is 8:1 in favor of the thin 
layer, which contributes to 89% of the precession ampli- 
tude (11% for the thick layer). These relative preces- 
sion amplitudes were expected from the relative weight 
of the thick and thin layers and from the approximate 
analytical model presented in section IIVBI The simu- 
lated field separation between the two coupled uniform 
modes {ojbqq -^Aoo)/i = 1-28 kOe compares also well 
with the 1.30 kOe estimate from the 2D model, with the 
dynamic dipolar coupling taken into account. Finally, 
one can check from the simulations the independence of 
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FIG. 11. (Color online) Simulated precession patterns of the 
eigen-vectors. Column (a) shows the precession profiles across 
the thin (light color) and thick (dark color) layers. Columns 
(b) and (c) show the dynamics in the thin Py^ and thick Pyt 
layers, respectively, with the color code defined in FIG. [6] 



the precession profiles on the thickness. This confirms 
the vaUdity of the 2D approximation and explains the 
performances of the analytical model. 

We now briefly comment on the simulated spectrum of 
FIG.[TOl which enables to identify the SW modes excited 
by the orthoradial Oersted field produced by the RF cur- 
rent flowing through the nano-pillar. From FIGS. [TUl and 
[TTl it is clear that the modes which couple to this exci- 
tation symmetry have a rotating phase in the azimuthal 
direction, characteristic of the i = +1 modes^l. We find 
that the SW modes of FIG. [TU] show the same series of 
A/B and n indices as those in FIG. IH] (but their i!-index is 
different). This sustains the translational correspondence 
between the SW spectra of FIGS. [2^ and[5}D. Finally, we 
point out that, for all the modes displayed in FIG. [TTl 
the pinning conditions at the boundaries of each disk are 
not trivial, which we attribute to the collective nature of 
the motion driven by the dipolar coupling2^. The general 
trend observed here is that the thin layer is less pinned 
than the thick layer for in-phase modes, and vice versa. 



To summarize, the 3D micromagnetic simulations en- 
able the identification (with three indices, A/B, i and 
n) of the SW modes probed experimentally by both a 
uniform RF magnetic field and an RF current flowing 
through the nano-pillar, i.e., of their respective selection 
rules. They confirm the experimental analysis performed 
in section IIIII and give a deeper insight on the collective 
nature of the magnetization dynamics in the nano-pillar 
discussed in section HVBl 



V. SYMMETRY BREAKING 

In the following, we review some characteristic spectral 
features associated with the breaking of the axial symme- 
try in our experiment. First, we experimentally report 
on the appearance of ^ = modes in the SW spectrum ex- 
cited by an RF current flowing through the nano-pillar, 
when a small tilt angle is introduced between the ap- 
plied field and the normal of the layers. This bridges the 
gap between our mechanical-FMR experiments and usual 
ST-FMR measurements^!^. Second, we have simulated 
the spectral distortions introduced by breaking the cylin- 
drical symmetry of the SW confinement potential. This 
enables to explain the lift of degeneracy in the SW spec- 
trum, which leads to the splitting of modes. 



A. Polar angle dependence 

The dependence on the polar angle 9h = (i,i/cxt) of 
the mechanical-FMR spectra excited by a uniform RF 
magnetic field and by an RF current fiowing through the 
nano-pillar is presented in FIG. [T^ Let us first focus 
on the conventional FMR spectra shown in FIG. [T^ . 
acquired at three different polar angles from the exact 
perpendicularity, increasing by steps of 1°. The main ef- 
fect here is the shift of the i = SW modes spectrum 
towards lower field as 0h increases, which has been ex- 
plained in details in Reff22, for a single magnetic disk. 
It is due to the decrease of the demagnetizing field pro- 
duced by the tilt of the equilibrium magnetization away 
from the normal. In fact, in each magnetic layer j = a, b, 
the uniform magnetization creates a non-uniform dipo- 
lar field AttMjG^'^'' *'"'ii which is maximum in the exact 
normal configuration. The equilibrium direction Uj is in 
the plane (i, i?ext) and makes a polar angle 9j > 9h with 
the normal determined by Eq. (|A34p . It can be estimated 
that when iJoxt ~ 10 kOe and 6h increases from 0° to 2°, 
the equilibrium angles 6a and 6h of the static magnetiza- 
tion in the thin and thick layers linearly increases from 
0° to Ri 9° and from 0° to « 13°, respectively. This leads 
to a shift to lower field of the FMR spectrum by about 
420 Oe (see appendix I A 2[) . in agreement with the data. 
We also emphasize that, in fact, the profiles of the SW 
eigen-modes are affected by the breaking of axial sym- 
metry, and that the pure £ = eigcn-modcs when 9h = 
become mixed with £ t modest when 9h + 0. 
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FIG. 12. (Color online) Dependence of the mechanical-FMR 
spectra excited by a uniform RF magnetic field (a) and by an 
RF current fiowing through the nano-pillar (b) on the polar 
angle 6h between the applied field and the normal to the 
layers. Superposed (in purple) is the behavior of the high 
field tail at larger power. 



We now turn to the influence of the polar angle 9h on 
the FMR spectra excited by an RF current {i^i = 170 /iA). 
The same global shift towards lower field as discussed 
above is observed in FIG. [12b by looking at the red spec- 
tra acquired with an increasing 6h- But there is an im- 
portant additional effect here. Whereas only £ = +1 SW 
modes arc excited by the RF current flowing through 
the nano-pillar in the exact perpendicular geometry, res- 
onance peaks can also be detected at the positions of £ = 
SW modes when 9h + 0. Although the amplitudes of the 
£ = modes are not large in FIG. Wlh . it is quite clear 
that they all grow as Oh increases. In order to reveal this 
effect better, we have reported in purple on the same fig- 
ure the resonance peak of the mode Bpo excited with a 
-\-\2 dB larger power [i^f ~ 680 /.tA), as a function of 6h- 



Despite the large RF current excitation, its amplitude al- 
most vanishes at Oh = 0. Then, it increases linearly with 
Oh ; until it becomes almost as large as when it is excited 
by the uniform RF field h,i ~ 2.1 Oe used in FIG. [Ha. 

The experimental data and their analysis presented in 
the previous sections [TT] to IIVI demonstrate that in the 
exact perpendicular configuration, only f. = +1 modes 
are excited by the RF current fiowing through the nano- 
pillar, due to the orthoradial symmetry of the induced 
RF Oersted field, Eq. (|A14p . Because there is no overlap 
between this particular excitation symmetry and the uni- 
form azimuthal symmetry of the £ = modes, the latter 
do not couple to the RF current excitation. The fact that 
these hidden modes in the exact perpendicular configura- 
tion can be excited by introducing a small misalignment 
angle between the applied field and the normal to the 
nano-pillar f is a striking result. It means that the selec- 
tion rules associated to the RF current excitation change 
if the applied field is tilted away from i, what we shall 
now explain. 

Due to the smaller demagnetizing field in the thin mag- 
netic disk than in the thick one (due to M^ < Af{,), the 
equilibrium angle of the thin layer is smaller than in the 
thick layer, Oa < 0^, as obtained from Eq. (|A34p . For 
the parameters of our nano-pillar, /3 = 0^ - Oa ~ 20 h, at 
Hcxt ~ 10 kOe and for a small angle Oh- It means that if 
Oh + 0, the magnetization vectors in both layers are mis- 
aligned from each other by an angle /3 = {Ma,Mb), so 
that the cross product Ua^Ub is finite and lies in the plane 
parallel to the layers, say along x. Thus, the spin transfer 
excitation (27rA)~^irtsin/3A associated to the RF current 
flowing through the spin-valve nano-pilla r^^'^^ , which is 
vanishing in the exact perpendicular configuration where 
/3 = 0, becomes finite if there is a small misalignment an- 
gle Oh t (see Eqs. (|A15[) and (|A16p in appendix lAli 
(27rA)"^ is the spin transfer efficiency). Because this so- 
called ST-FMR excitation has the same symmetry as an 
in-plane uniform RF magnetic field, it is expected to ex- 
cite SW modes having the £ = 0-index symmetry. Still, 
this excitation has to compete with the RF Oersted field 
excitation, which is independent of Oh and is much larger 
in our configuration due to the small value of /3 (< 5°). 
Therefore the amplitudes of the £ = +1 modes are much 
larger than those of the £= modes in FIG. [T^ . 

It is also clear that the amplitude of the mode i?oo ex- 
cited by the RF current (purple peaks in FIG.[T^) grows 
linearly with Oh, as expected from the above discussion. 
We emphasize that a quantitative understanding of the 
amplitude of the peaks excited by ST-FMR would re- 
quire to consider the collective nature of the dynamics 
in the nano-pillar and the asymmetry of spin transfer in 
the thick and thin magnetic layers. Finally, we note that 
the small signal observed at Oh = should in principle 
vanish with the ST-FMR excitation. This reminiscent 
signal can be ascribed to a small misalignment of the ap- 
plied field with respect to the normal to the nano-pillar 
(the precision on the orientation is 0.2°) or to a slight 
asymmetry of the RF current lines through the nano- 
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FIG. 13. (Color online) Simulated SW spectra for a nano- 
pillar with an elliptical section (see text). Linear response to 
a homogeneous RF magnetic field excitation (a) and to an 
orthoradial RF Oersted field excitation (b). The precession 
patterns of the lowest energy modes are shown in the insets. 



pillar, which would induce a small asymmetry of the RF 
Oersted field, thereof adding a small in-plane uniform 
component to the orthoradial magnetic field. 

To summarize, this study enables to derive the selec- 
tion rules of the RF current excitation. In the exact 
perpendicular configuration, the magnetizations of both 
layers are aligned, and only i = +1 modes can be excited 
due to the orthoradial symmetry of the current-created 
Oersted field {£= modes are hidden) . But when a finite 
angle is introduced between the magnetizations in each 
layer by slightly tilting the applied field away from the 
normal, £ = modes can be excited by ST-FMR, which 
has the same symmetry as a uniform RF field excitation 
polarized in-plane. 



B. Confinement asymmetries 



As seen in section llVCi the 3D micromagnetic sim- 
ulations enable to identify the SW modes observed in 
the experimental spectra. Still, the latter are more rich 
than the simulated power spectra, due to the splitting of 
some resonance peaks, which was noted in sections III CI 
and nil M In particular, the experimental peak at Ho, 
identified as the mode -Bio, is clearly split in two, with 
a smaller resonance about 100 Oe away in the low field 
wing of the main peak, which is not the case in the sim- 



ulation (see FIG. [10)) . The peaks at H@ (identified as 
Bii) and at Hq (Aiq) are also split, contrary to the sim- 
ulations, where all these peaks are single. In contrast 
to these observations, the peak at ff®, which is the uni- 
form mode Bqq, is single both in the experimental and 
simulated spectra. 

So, it seems that experimentally, the occurrence of the 
mode splitting depends on the mode index, whereas in 
3D simulations, in which the nano-pillar has a perfect 
cylindrical shape, none of the resonance peaks is split. 
This suggests that the observed splittings are related to 
asymmetries in the confinement of the disks, and that the 
various SW modes are affected differently because they 
probe different regions. The fact that the double peak at 
Ho depends on the tilt angle (see FIG. \T^) and is more 
or less pronounced depending on the direction in which 
the applied field is tilted from the normal (not shown) is 
another strong indication that some symmetry breaking 
in the lateral confinement is at the origin of this efFecti2^. 

To support this idea, we have carried out new 3D mi- 
cromagnetic simulations with the SpinFlow 3D package 
on a structure that break the perfect cylindrical sym- 
metry of the nano-pillar. We have kept a perfectly flat 
structure, but we have used an elliptical cross-section. 
The long axis of the ellipse is 250 nm, while the short axis 
is 200 nm. The influence of this breaking of symmetry on 
respectively the i = spectrum (RF field excitation) and 
the £ = +1 spectrum (RF current excitation) is presented 
in FIG. m 

Concentrating first on the standard FMR SW spec- 
trum of FIG. [T5k . one can see that the lowest energy 
mode -Boo remains a single peak. This illustrates the 
intuitive idea that the uniform SW mode, where the os- 
cillation power is mostly concentrated at the center, is 
rather not sensitive to change of the confinement at the 
periphery. The same behavior applies for the lowest en- 
ergy mode of the thin layer, Aqq. 

The simulated SW spectrum of FIG. [T^b shows a dif- 
ferent sensitivity to the shape asymmetry. It is observed 
that the lowest energy mode with the £ = +l-index splits 
in two peaks, in contrast with the lowest energy £ = 
mode which remains a single peak. We also note that 
the satellite peak, induced by the elliptical confinement, 
is located in the low field wing of the main resonance, as 
in the experiments of FIG. [^b- The precession patterns 
shown in the two insets FIG. [T5b reveal that the ellip- 
tical shape introduces some mixing between the £ = +1 
and £ = -1 SW modes (the £ = -1 mode corresponds to 
SWs that are rotating around the disk in the opposite 
direction as the Larmor precession). In a circular disk, 
these two modes are degenerate, and only the £ = +1 
mode couple to the orthoradial Oersted field excitation. 
But in the ellipse, the two eigen-modes split and become 
mixed, as shown by the two eigen- vectors displayed in 
FIG. [T3b . which correspond to the linear combinations 
Ji{p) + J-i{p) a-nd Ji(/o) - J-i(/o). The simulated spec- 
trum of FIG. 113b reproduces well the main features of 
the mechanical-FMR spectrum of FIG.[5]3, including the 



18 



splittings observed for the £ = +1 SW modes (revealed 
even better by injecting a do current through the nano- 
pillar, see FIG. [3}d)- Thus, a small elliptical shape pro- 
duced during the nanostructuration of the nano-pillar is 
most likely responsible for the double peak observed at 
Ho in FIG. Eb- 

To summarize, the comparison between 3D simulations 
and experiments demonstrate that the observed mode 
splittings originates from a small asymmetry in the lat- 
eral confinement of the nano-pillar. 



VI. CONCLUSION 

In summary, we used the MRFM technique^ to study 
the SW eigen-modes in the prototype of a STNO - a nor- 
mally magnetized nano-pillar composed of two magnetic 
layers coupled by dipolar interaction. 

In contrast to transport spectroscopy techniques^—, 
MRFM is sensitive to all SW modes excited in the 
sample^ and is completely independent of the transport 
properties of the studied spin-valve sample. Therefore, 
MRFM provides an alternative and complementary view 
on the magnetization dynamics in hybrid magnetic nano- 
structures. The additional advantages of the MRFM 
technique are its high sensitivity (in this study, it was 
able to detect angles of precession as low as 1° in the thin 
magnetic layer) and its ability to operate on standard 
STNO devices buried under contact electrodes without a 
specific probe access to the studied sample^. 

Using MRFM, we were able to compare the SW spec- 
tra of a passive perpendicularly magnetized STNO-like 
sample excited by a uniform in-plane RF magnetic field 
and by an RF current flowing perpendicularly through 
the layers. We found that distinctly different SW modes 
(having azimuthal indices £ = and £ = +1, respectively) 
are excited by the two above mentioned excitation meth- 
ods. By studying the influence of a spin polarized dc 
current on the observed SW spectra we were able to 
determine which of the magnetic layers of the studied 
nanopillar plays the dominant role in the magnetization 
dynamics resulting in the appearance of each particular 
SW mode. 

We also developed a simple analytic theory allowing 
to perform a comprehensive labeling of all the SW eigen- 
modes of a magnetic nanopillar in the studied axially 
symmetric case. This labeling requires three independent 
indices: the usual azimuthal and radial indices £ and n 
used for the SW modes of a single magnetic disk and an 
additional index referring to the binding or anti-binding 
(B or A) coupling between the two magnetic disks form- 
ing a nanopillar. The obtained experimental and ana- 
lytic results were also compared to the results of 3D nii- 
cromagnetic simulations obtained with the SpinFlow 3D 
package^, which confirmed the mode labeling obtained 
from the analytic theory. 

Thus, we learned that in the axially symmetric case 
of a perpendicularly magnetized nano-pillar, the excita- 



tions by the RF field the RF current lead to two mutually 
orthogonal (and mutually exclusive) sets of excited SW 
modes: only the £ = modes are excited by the uniform 
RF magnetic field, while only the £ = +1 modes are ex- 
cited by the RF current. Therefore, the ^-index, related 
to the azimuthal symmetry of the SW modes, is the dis- 
criminating parameter for the selection rules of the SW 
mode excitation. 

Moreover, we have demonstrated experimentally and 
numerically that the mode selection rules are affected 
by the breaking of the axial symmetry of the studied 
nano-pillar, either by tilting the bias magnetic field or by 
making the sample cross-section elliptical. In particular, 
if the axial symmetry is broken by tilting the bias mag- 
netic field, the £= modes can also be excited by an RF 
current. This excitation is caused by the ST-FMR mech- 
anism working when the magnetization vectors in the 
two magnetic layers of the nano-pillar arc not coUinear. 
Also, the importance of the dynamic dipolar interaction 
between the magnetic layers of the nano-pillar have been 
clearly demonstrated by our results. 

We believe that our results are important for the opti- 
mization of the characteristics of nano-spintronic devices, 
and in particular STNOs, and for the experimental de- 
termination of the STNO parameters. 

First of all, an accurate identification of the SW modes 
that can be excited in an STNO nano-pillar is necessary 
to understand the details of the high frequency STNO 
dynamics. The proposed identification of the nano-pillar 
SW modes can be used for the experimental determi- 
nation of the nano-pillar characteristics, such as radius, 
static magnetization, gyromagnetic ratio, and dissipa- 
tion. We note that in traditional STNO experiments, 
where the magnetization dynamics in a magnetic nano- 
pillar is excited by a spin-polarized bias current creating 
a significant Oersted magnetic field with the £ = +1 sym- 
metry, it is easy to mix-up the £ = and the £ = +1 SW 
eigen-modes. The spectra of these modes arc in almost 
translational correspondence, and the experimentally ob- 
served dependence of the mode frequencies on the bias 
magnetic field can be well described by the traditional 
Kittel expression (see e.gM-). Thus, the possibly excited 
£ = +1 mode can be easily interpreted as a ^ = mode, 
which will lead to the apparent reduction of the "free" 
layer static magnetization necessary to fit the Kittel ex- 
pression for the mode frequency. For instance, mislabel- 
ing the lowest energy mode of the £ = +1 SW spectrum as 
the uniform mode (€, n = 0, 0) combined with a small mis- 
alignment of the applied field would lead to a discrepancy 
as large as 1 kG in our case. 

Second, the fact that in most cases both magnetic lay- 
ers of a nano-pillar take part in current-induced mag- 
netization dynamics is very important for the correct 
identification of the excited SW modes. The collective 
(coupled) character of the SW modes in a nano-pillar 
can directly influence the magnitude of the spin trans- 
fer torque, which is dependent on the relative orientation 
of the magnetization vectors in the two magnetic lay- 



19 



ers. One might expect, that the efSciency of the spin 
transfer torque for a particular SW mode depends not 
only on which layer ( "free" or "fixed" ) is dominating the 
mode dynamics, but also on the coupling (in phase or 
anti-phase) between the magnetization precession in two 
layers. In our experimental case, the interlayer coupling 
is in-phasc for the SW modes dominated by the dynamics 
in a "free" (thin) layer. To obtain an interlayer coupling 
which is anti-phase for the SW modes dominated by the 
dynamics of the thin layer would require for example to 
increase its magnetization compared to the "fixed" layer 
one. 

Finally, it is important to note that the MRFM tech- 
nique has allowed us to study spin transfer effects in 
the axially symmetric configuration of a perpendicularly 
magnetized nano-pillar, where the excitation of magne- 
tization dynamics by ST-FMR vanishes due to the sym- 
metry reasons. This geometry is rather important for 
applications as the excited SW modes have the maxi- 
mum non-linear frequency shift coefficient^. This cre- 
ates the maximum agility of the mode frequency with 
the bias current and, therefore, the maximum width of 
the synchronization band to the external periodic signal 
and to the large arrays of other STNOa^^. Phase syn- 
chronization has been identified as a possible mean to 
dramatically increase the generated microwave power of 
these nano-oscillators and, at the same time, reduce their 
HnewidthSS,-^. 
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Appendix A: Theoretical material 

1. Equation of motion 

In this appendix, we detail the derivation of the equa- 
tions which govern the dynamics of a ferromagnetic layer 
in the presence of an external periodic excitation and 
of spin transfer, following the general formalism intro- 
duced in section IIVAII For an isolated layer, the local 
dynamics (within the exchange length) of the magneti- 
zation vector is described by the Landau-Lifshitz (LL) 
equation: 



IdM 

7^r 



HxM + h{t)xM, 



(Al) 



where the second term on the right-hand-side of Eq. (JA1[) 
represents the perturbation term. The field H is the ef- 
fective magnetic field: 



H= Ho- 47rG * M 



(A2) 



Here, Hq is the total static external magnetic field (possi- 
bly spatially-dependent) and the linear tensor self-adjoint 
operator G describes the magnetic self-interactions. 

Considering only the linear processes, we can repre- 
sent the time dependent (out-of-equilibrium) part of the 
magnetization as a series over the SW eigen-modes: 

M(t,r)-Msu{r)«m{t,r) (A3) 

~ Y.c,.{t)m,,{r) +C.C. , 

where Ms is the saturation magnetization of the layer. 
Here, c.c. stands for the complex-conjugated part. The 
coefficients c^it) are time-dependent SW amplitudes. 

The second term on the left-hand-side of Eq. (|A1[) 
represents the perturbations from the equilibrium state, 
including the non-adiabatic contributions. The non- 
conservative perturbation magnetic field h(fy may de- 
pend on time and be a function of the magnetization 
distribution M(t). It can be approximately represented 
as: 



h{t) = hi{t) +Li *m{t), 



(A4) 



where hi(t) is the external perturbation field and Li is 
a certain linear operator, allowing h{t) to depend on the 
magnetization distribution. The latter case may describe 
the influence of the Gilbert damping -{a/jMs)dM/dt 
through: 

Li*m = i— y^ u};,{c,y{t)m^{r) - c.c.) , (A5) 

7 u 

where a is the Gilbert damping constant. 

Substituting the series representation Eq. (|A3|) and the 
representation of the perturbation field Eq. (jA4[) into 
Eq. (jAip and using orthogonality relations Eq. ([7]), one 
can obtain the following equations for the SW amplitudes 

— ^ = -iuj„c^ + ijY.(S^^^'C„' + S^yCi,,\ + i^hy , (A6) 

at j^r 

where 



with 7 being the modulus of the gyromagnetic ratio. 
The LL equation is written here in its perturbative form, 



^ly.i^' 



S..-7- 



hu 



[m^ 


■Li *m^,) - {{u-hi){m^-m^,)) 




N, 


{rriu 


■Li *rn^f) - {{u-hi){rn^-rn^,)) 



N. 



{my -hi) 



(A7a) 

7 

(A7b) 
(A7c) 
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In many cases the perturbed equations (|A6p can be 
further simpUfied by retaining only the diagonal term 
SvM = iauj^{rfi^ -m^) /{jJV,y) (assuming that there are 
no degenerate modes). The SW damping rate is then 
given by 



r^ = acj. 



AC 



(A8) 



The damping rate r^^ is responsible for the finite 
linewidth of the resonance peaks, AiJ (FWHM). If the 
sample is homogeneously magnetized and the precession 
is circular, the simple relation aAH = uj^h holds. 

From the equations above, one recovers for the coeffi- 
cient c^ the equation of motion of a damped harmonic 
osciUator: 



dt 



-iw,yCy - T^Cy + i^hy . 



(A9) 



If a second magnetic layer j' is electrically connected 
to layer j and spin transfer is allowed between them, the 
equation of motion must be modified. When a charge 
current I is flowing through the layers, the additional 
Slonczewski-Berger term^ii adds to the operator Li of 
Eq. (|A4[) . It is possible to define an effective damping in 
layer j, which depends on the spin polarized current and 
on the relative angle between the magnetization in layer 
j and the direction of the spin polarization Uji^: 

^ {fn„-my) le {{rfiu ■mu){uj ■ Uj')) 



where 



2eN, 



N,- ^ ^ 



M, 



•yh 



(AlO) 
(All) 



is the dimensionless total number of magnons that can 
be excited inside the volume Vj of layer j. Here, h is the 
reduced Planck constant, e the modulus of the electron 
charge and e the spin polarization efficiency of the cur- 
rent. The threshold current for auto-oscillations in layer j 
corresponds to F* = 0, i.e., if Uj j| u'-, I^h = -^aUyNgeje. 
Using e = 0.3 and the parameters of our thin layer, one 
can estimate /th - -4.8 mA for the uniform SW mode 
at 8.1 GHz, in agreement with the experimental datai^^. 
This result is also in quantitative agreement with calcula- 
tions performed for our nano-pillar device with no adjust- 
ing parameters in the framework of continuous random 
matrix theory (CRMT) described in Ref.^^. 

We now turn to the periodic external excitation hi (t) = 
hi exp*'^* , whose amplitude hi is composed of three dif- 
ferent contributions, 



hi = h^ + hoc + h 



ST : 



(A12) 



that we shall detail below. 

The first type of excitation corresponds to a uniform 
RF magnetic field applied perpendicularly to the effective 



field H. This configuration corresponds to conventional 
FMR spectroscopy. Assuming that the sample is uni- 
formly magnetized along the nano-pillar symmetry axis 
f , it reduces to: 



h^ = h^fX, 



(A13) 



where A is a unit vector in the in-plane direction and /i,f 
the linearly polarized amplitude. 

SW spectroscopy can also be performed by injecting 
a uniform RF charge current ij-f through the nano-pillar 
(i.e., along z). First, this produces an orthoradial RF 
Oersted field: 



h 



Oc 



47r 

To 



_hi_p_ 

2ttRR 



{-siiKpx + cos0y) 



(A14) 



where R is the radius of the nano-pillar and (p, (^) are 
the polar coordinates. In this formula, the current should 
be expressed in A and the prefactor between the square 
brackets converts A/cm into Oe (cgs units). The maxi- 
mum amplitude of the RF Oersted field is reached at the 
periphery of the nano-pillar, p = R, and equals 1.6 Oe 
for a peak amplitude irf = 100 fiA and the experimental 
parameters. 

Second, the RF current produces a ST-FMR excita- 
tion: 



hsT 



2ttX 



'] 



(A15) 



where we have rewritten the spin-transfer efficiency of 
the charge current in Eq. (|A10[) as a function of 



27rA = 7 



2eiV,, 



(A16) 



which has the dimension of a distance (A ^ 200 nm for 
our thin layer). If the thin and thick layers arc mis- 
aligned by an angle /3 in the plane {y.,z), Eq. (jA15|) re- 
duces to ft-sT = irt/(27rA) sin/3A, which demonstrates that 
the ST-FMR excitation is equivalent to a linearly polar- 
ized RF magnetic field, Eq. (|Il3)) . The ST-FMR exci- 
tation vanishes if the magnetic layers are parallel. The 
amplitude ratio between the ST-FMR and the RF Oer- 
sted field excitations, both produced by the RF current 
flowing through the nano-pillar, is /isx/ZiOe - (-R/A) sin/3. 
In our geometry, A r* 2i?, but due to the small angle (3 
between the layers, the RF Oersted field contribution is 
much larger than the ST-FMR one. We note that even if 
j3 M 7r/2, the contribution of the RF Oersted field cannot 
be disregarded in general in ST-FMR experiments. 



2. Numerical application 

In this section, we derive a practical guideline to cal- 
culate the eigen-frequencies w using the analytical for- 
malism developed in section HV Al Let m^ be a certain 
orthogonal basis in the space of the vector functions m 
satisfying both the local orthogonality to ii and the total 
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pinning condition at the boundary of the magnetic body. 
Spin-wave eigen modes can thus be expressed as a series 
expansion on the rUj^ basis (cf. Eq. ()A3|) ). A general ex- 
pression for the eigen-frequencies can be found from the 
condition of vanishing determinant: 



\N^'^^u!-Afi,',i,{fi}^, ^ II = 



(A17) 



where for the simphcity of the discussion, it is convenient 
to introduce a curly bracket notation, to indicate that 
the enclosed quantity is spatially weighted by the spatial 
pattern of the mode profile and averaged: 









(A18) 



This echoes the chevron bracket notation introduced in 
Eq. ([1]) to indicate the homogeneous spatial average over 
the volume of the magnetic body. Here Af represents a 
renormalization quantity, defined by 



Ny'^i, = i{my, ■ {u X m^)) 



(A19) 



which has in general off-diagonal elements. 

In the case of perpendicularly magnetized disks, where 
the set of Bcssel functions ^{x + iy)e^^^'^ j£{k£,np) diago- 
nalizes the uncoupled Hamiltonian, the secular Eq. (|A17[) 
becomes diagonal and we recover Eq. ([8]): 



m. 



(A20) 



We shall now perform the numerical application of the 
eigen-value of the lowest energy mode {i,n = 0,0) using 
the parameters of our nano-pillar shown in Table HI We 
will drop the subscript ;/ to the curly brackets, under- 
standing that the spatial average in Eq. (|A20p is made 
over the uniform mode rrii, = ■^Jo{kop){x + iy), where 
fco = 2.4048/i? is its wave-vector. In this case the value 
of the normalization constant is simply A/q = (Jo) ~ 
Ji(koR) = 0.2695. 

The different contributions that enter inside the oper- 
ator {rj} are detailed in Eq. (|3|): 



{n} = -f{H}+47r-fMj{G} . 



(A21) 



We start with the calculation of the amplitude of effective 
magnetic field, the first term on the right hand side of 
Eq. (|A2ip . As shown by Eq. ([5]), the scalar value H 
along z decomposes itself in two terms: 



{H} = {z ■ Ho} - 47rA/^ {z ■ G * z } . 



(A22) 



The term z ■ G * z represents the static magnetic self- 
interaction. In the case of homogeneously magnetized 
body, the inhomogeneous exchange contribution to the 
static self-interaction is strictly null and the second term 
of Eq. (|A22[) reduces to the magneto-dipolar contribution 
G^'^\ which has the following form in the wave- vector 
representation: 



G('^)(r) = |i?(fc)^exp*'^ci3fc, 



(A23) 



where D{k) is the Fourier transform of the body shape 
function's and the symbol ® denotes direct product of 
vectors. For a disk of radius R and thickness i, an an- 
alytical expression for the different position-dependent 
demagnetization tensor elements of a disk Nu^[ii^fi{r) = 

u{r) ■ G^"^' * V valid in the whole space are available 
in Ref.— . For perpendicularly magnetized disks where 
u = i) = z, the expression of the self- integral becomes 



{^i"^} = 7727 [dU^ Joikop)N,z[R,t,]{p,z) : 



(A24) 



for both the thin {j = a) and thick (j = b) layers. Their 
numerical values arc displayed in Table. IIIII 

The term {z ■ Hq) of Eq. (jA22|) is the projection on 
the precession axis of the total applied magnetic field. 
It comprises the external magnetic field {-ffoxt} = ^ext, 
the stray field of the mechanical-FMR probe {i?sph} = 
190 Oe and the cross-magneto-dipolar static interactions 
between each layer. The latter can be estimated from the 
cross tensor elements of the static magneto-dipolar field 
of the j'-th disk produced over the volume of the j'-th 
disks: 

{^-''''} = 772Y f ^'^^^ Jl{kop)N,, [R,t^,]{p, z + zo), 
Wo / ^3 

(A25) 
where zq is the distance between the centers of the two 
axially aligned disks. The numerical values of the cross 
tensor elements are reported in Table. IIIII Putting all the 
above elements together, the total effective field simply 
writes: 

{H} = iJext + {i^sph} - 47r {TVi^'^)} M, - An {n^^^'^} M,, . 

(A26) 
We now turn our attention to the integration 
AttMj {G}, the second term on the right hand side of 
Eq. (|MT|) . We recall that for Py the operator G = G'^^^ + 
G^'*' is the sum of the inhomogeneous exchange and 
magneto-dipolar interactions. In the wave-vector rep- 
resentation, G^*^-* = Al^k^I, where the exchange length 

Acx = \/2^/(47rA/^) depends on the exchange stiffness 

constant J, expressed in erg/cm (= 10"^ in Py). It pro- 
duces the exchange field: 



{i^ox} 



AnMjAl^k^ , 



(A27) 



which yields the value {Hex} = HO Oe. 

The other contribution is the dynamic magneto-dipolar 
self-interaction, which represents the depolarization field 
of the SW mode on itself. For the i = modes, an ana- 
lytical expression can be derived: 



{""-'} ^Th)I/'''''^'''^L 



R dN^ 
du 



:[u,tj, 



]{P^z) 



du 



Joikou) 



(A28) 

where the quantity in the second integral is the magnetic 
stray field produced at the spatial position r by a cylin- 
drical tube of width du, radius u, and thickness tji, ho- 
mogeneously magnetized^^ along x by Jo{kou). We use 
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the same expression above to write the self- and cross- 
contribution, understanding implicitly that the spacer 
value zq should be added in the later case, as shown in 
Eq. (|A25p . The values of the self- and cross-tensor ele- 
ments are reported in the last line of Table. IIIIl We men- 
tion, that an approximate expression of the self-Eq. ()A28|) 
has been derived by Kalinikos and Slaving for the lowest 
SW branch of platelet shape bodies with uniform mag- 
netization across the film thickness. This expression re- 
duces to 



{N..}-^{l-Gt) 



(A29) 



where the analytical expression of G^ „ for Bessel func- 
tions is given by Eq. (26) in re^. 

The cross elements are responsible for the dynamic 
dipolar coupling detailed in section IIVBI 



b'j' 



47r 



{n!^;''^} 



M. 



j' 



(A30) 



The value of the coupling frequency fl for the lowest en- 
ergy mode yields: 



Q' 



7 



{Nf^ 



} 4^A/fc JTvi^'"^ } AnMa , (A31) 



which leads to J7/27r ~ 0.56 GHz. 

Neglecting the dynamical dipolar coupling (the gener- 
alization to 17 + is Eq. (jisp). we derive an expression 
for the eigen- value of index v = jo,o- 



7 



{iJ}+47r{iV(^/")}A/,+{iJex}, (A32) 



j = a,b being the layer index and {H} being defined in 
Eq. (jA26[) . Eq. (|A32p is a simplified expression valid for 
circularly polarized modes (£ = index), where we have 
taken advantage of the equality {Nj:x}o - {Nyy}o in our 
circular disk. This expression can be extended to higher 
orde r mod es by using {Nxx}e,n ~ {Nxx}o\/ke.n/ko inside 
Eq. (jA32p . This approximation is derived from the el- 
lipticity of £ + modes {{Nxx} + {^yy})- One needs 
thus two separate equations (|A32[) for the values of uj for 
each cartesian axis^: one proportional to rn^Knixiny), 
the other to my/{mxmy). The product of these two 
equations is independent of the ellipticity, leading to the 
general expression for the eigen- value of arbitrary index 
'^ = je,n- 

{{H}^ + 4^{^^^'^^}, M, + {^cx} J , (A33) 

which can be seen as a generalization of the Kittel for- 
mula for arbitrary shaped multi-body. 

Equating Hagg = H® and Hbgg = H® in Eq. (|A32p . 
where Hagg and H^gg are the resonance fields at /fix = 
8.1 GHz of the uniform modes in the thin and thick disks, 
respectively, leads to AnAda = 8.0 x 10'^ G and AttMi, = 
9.6x103 G. 



TABLE III. Values of the self- and cross- depolarization tensor 
elements weighted by the precession profile of the uniform 
mode for the thin (j = a) and thick (j = b) disks. 





(a, a) 


(a, 6) 


ib,a) 


(6,6) 


{iVii'' ^} 


+0.979 


-0.068 


-0.017 


+0.919 


{iv£-^'^} 


+0.016 


+0.042 


+0.011 


+0.056 



Finally, the above formalism also enables to determine 
the angle 9j = {z,Uj) between the equilibrium direction 
of the magnetization in layer j and the normal axis when 
the bias field is applied at a polar angle 9h = (i,H'cxt)- 
From the equality u x H^xt = 47r7\/jM x I^G^'^'' >i- u), one 
extracts the relationship: 

i/ext Sin (% - Oh) = 27tMj {{Nii^^) - (N^i'^)) sin 20, 

(A34) 
The angle 9j is useful to estimate the shift to lower 
field of the FMR spectrum^, 2TTMj{{Nii>}-{N^i^}){l- 
cos26'j) « 420 Oe, when Ot = 13° in the thick layer. 



Appendix B: Methods and calibration 
1. Mechanical vibration amplitude 

Here, we detail the experimental protocol used to cal- 
ibrate the amplitudes of the uniform RF magnetic field 
and of the mechanical-FMR signal. The procedure uses 
the non-linear properties of the magnetization dynam- 
ics and consists in studying the power dependence of the 
line shape. In the following, AH denotes the FWHM 
linewidth measured in the linear regime. 

We use the onset of foldover as a mean to calibrate the 
strength of the RF field produced by the microwave an- 
tenna. This non-linear effect is responsible for the asym- 
metric shape of the resonance peaks in FIGS. [5] and |31 
In fact, it was pointed out by Anderson and Suhl^ that 
the resonance curve at high power should be skewed, due 
to the static change of the magnetization AI^, which also 
shifts the resonance frequency. For a normally magne- 
tized sample, this non-linear frequency shift is positive 
(blue-shift), and the field-sweep line shapes are distorted 
towards low field. There is a critical strength of the RF 
magnetic field he (linearly polarized amplitude) for which 
the slope of the resonance curve becomes infinite on the 
low field side of the resonance^: 



2AHy 



2AH 



\ 3V3\{N,,} - {Nxx}\^7rM, 



(Bl) 



where {N^z} - {N^x} is the difference between the de- 
polarization factors in the longitudinal and transverse 
directions. Experimentally we find that for the peak 
at Hq), this onset is reached when the output power of 
the synthesizer at 8.1 GHz is Pq = +9 dBm. Using the 
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magnetic properties of the thick layer (Table |I| , we infer 
from Eq. (|Bip that at the critical onset of foldover, the 
strength of the RF magnetic field is he = 4.2 ±0.8 Oc. We 
note, that this value is in agreement with the estimation 
made by directly evaluating the field produced by the RF 
current flowing in the antenna at 8.1 GHz for this output 
power, hj-f = 5.5 ± 1 Oe. 

Furthermore, this procedure gives a calibration of the 
amplitude of the mechanical-FMR signal (AM^). At the 
onset of foldover, the longitudinal change of the magne- 
tization is indeed^ 



47r(AA/,) 



3\/3 



AH. 



(B2) 



A numerical application of Eq. (jB2p yields 4:Tt{AMz) = 
36 ± 4 G, which corresponds to the critical angle of pre- 
cession {9c) = 5°. We have used this calibration of the 
cantilever vibration amplitude to evaluate the change of 
the longitudinal magnetization at the maximum of the 
peak at if® in FIG. Hi. 



2. Microwave setup 

In this appendix, we give some details on the mi- 
crowave circuit, which was carefully design to minimize 
the cross-talk between the RF field and RF current exci- 
tation parts. 

The calibration of the RF magnetic field produced by 
the microwave antenna has been presented in the previ- 
ous appendix. In order to calibrate the RF current flow- 
ing through the nano-pillar with respect to the synthe- 
sizer output power injected into the contact electrodes, 
we have first used a standard microwave setup. The 
nano-pillar electrodes are directly connected to the mi- 
crowave synthesizer through a picoprobe, a bias-T, and 
a semi-rigid coaxial line, that allow to perform voltage- 
FMR spectroscopy. In this experiment, the amplitude of 
irf flowing through the nano-pillar can be accurately de- 
termined, owing to the determination of losses and reflex- 
ions in the microwave circuit using a network analyzer. 
Then, the same experiment is repeated inside the MRFM 
setup, in which the contact electrodes are wire bounded 
to a microwave cable and the circuit contains more con- 
nections. The comparison with the standard setup yields 
an estimation of the rms amplitude of the RF current 
in the mechanical-FMR setup: i^f = 170 ± 40 /zA for an 
output power of -22 dBm injected at 8.1 GHz through 
the contact electrodes. 

It is also possible to estimate experimentally the high 
frequency coupling between the microwave antenna and 
the electrodes that contact the nano-pillar. For this, we 
exploit the fact that in the exact perpendicular configura- 
tion, different SW modes are excited by the uniform RF 
field {£ = 0- index) and by the RF current {£ = +l-index). 
If the RF magnetic field used to excite the conventional 
FMR spectrum would induce any relevant RF eddy cur- 
rent through the nano-pillar, £ = +1 modes which are 



excited in the SW spectrum of FIG. [2]3 should also be 
detected in the SW spectrum of FIG. [5^, which is not 
the case. We deduce from this observation that for an 
output power of +3 dBm injected in the antenna, the in- 
duced eddy current through the nano-pillar is less than 
when injecting -38 dBm directly through the contact 
electrodes, i.e., ?rf < 30 ^A. So, at /gx = 8.1 GHz, the 
isolation between the two parts of the microwave circuit 
is better than 40 dB. However, we note that the latter 
depends on the frequency, and that for some particular 
values, it can drop to only 20 dB. 

Still, owing to the broadband design of the contact 
electrodes and to the low microwave power required to ex- 
cite SW modes with the RF current excitation part, it is 
possible to acquire FMR spectra at a fixed bias magnetic 
field i/fix by sweeping the frequency of the RF current. 
We mention that in the frequency-sweep experiments pre- 
sented in FIG. [3^, the output power of the synthesizer 
is kept at -22 dBm over the full frequency range (4 to 
18 GHz), which results in an amplitude variation of i^f, 
mainly associated to frequency dependent losses in the 
circuit. We also note that the same frequency-sweep ex- 
periment cannot be performed as cleanly with the RF 
field excitation due to the high power that has to be in- 
jected in the microwave antenna and to the dependence 
of the isolation on frequency mentioned above. 



Cavity-FMR characterization of the extended 
film 



Before the nano-fabrication of the nano-pillar devices, 
a reference film of Cu60 | PyblS | CulO | Pya4 | Au25 
(thicknesses in nm) is cut out from the Si wafer for char- 
acterization purpose. The extraction of the material pa- 
rameters is obtained independently on this reference film 
by a reflexion X-band spectrometer (9.6 GHz) operat- 
ing at room temperature. The experiment consists in 
measuring the resonance spectra of the multi-layer as 
a function of the polar angle 9h between the applied 
field and the normal to the film. The resonance field of 
the layer j = a,b as a function of 6h depends only on 
the gyromagnetic ratio 7 and on the total perpendicular 
anisotropy field, which here reflects entirely the demag- 
netizing field A-kMj of the layer—. The obtained values 
for the gyromagnetic ratio (identical for both layers) and 
the magnetizations arc collected in Table ID The magne- 
tization of the thick layer (9.6 kOe) corresponds to the 
expected value for bulk Py with composition Ni8oFe2o- 
The magnetization of the thin layer is 1.4 kG smaller, 
which refiects the reduction of the magnetization in the 
interfacial layer (of the order of 1 nm) , due to the gradual 
composition variation of the NiFe alloy from Ni8oFe2o to 
the normal metal (Cu or Au)^*^. 

An estimate of the damping parameter and the amount 
of inhomogeneous broadening can also be obtained from 
the angular dependence of the linewidth of the resonant 
mode associated to each layer. The linewidth is in gen- 
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eral the sum of two contributions: an intrinsic relaxation 
of the magnetization vector (homogeneous width) and 
an inhomogeneous broadening corresponding to a dis- 
tribution of resonance fields (whose main sources have 
been described for a Permalloy polycrystalline layer—). 
The intrinsic damping parameter is deduced from the 
parallel geometry linewidth. In the Py;, 15 nm thick 
layer, the line widths observed in the parallel and per- 
pendicular geometries are respectively: A7f||(, = 64 Oe 
and AHj^iy = 73 Oe. The higher value observed in 
the perpendicular geometry reveals a fair amount of 
inhomogeneities^^. The linewidth observed in the paral- 
lel geometry corresponds to an intrinsic damping param- 
eter ab = (0.9 ± 0.1) X 10"^. For the thin Pya 4 nm layer, 



AH^a = 83 Oe and AHj_a = 171 Oe. From the parallel 
geometry linewidth, we deduce an intrinsic damping pa- 
rameter aa = (1.5 ± 0.3) X 10"^. This value is higher than 
for the thick layer because of the larger effect of the dif- 
fusion of the microwave magnetization of the conduction 
electrons in the adjacent normal metal layers^, associ- 
ated to the fact that the thin layer thickness (4 nm) is 
less than the spin-diffusion length in Py. The much larger 
value AH la is associated to a large contribution of the 
inhomogeneous broadening arising from a substantial ef- 
fect of intcrfacial roughness (and dispersion of heights of 
the crystallites of the base) on the thin Py layer grown 
on top of a 85 nm thick metallic base. 
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